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Preface 


The  goal  of  this  study  was  to  provide  a  single  source  document  on  the  effects  of  noise  jamming 
on  pulsed  radar  receivers.  This  study  would  cover  three  main  topics  in  depth.  These  areas  to  be 
covered  are  the  sources  of  electrical  noise,  various  noise  jamming  waveforms  and  the  effects  of  each 
jamming  waveform  on  various  pulsed  radar  receivers. 

This  thesis  effort  is  concentrated  on  determining  the  power  spectra  of  various  noise  jamming 
waveforms,  determining  the  pr  lability  densities  of  the  outputs  of  conventional  pulsed  radar  re¬ 
ceivers  and  the  effects  of  CW  interference  on  conventional  pulsed  receivers.  As  such,  this  thesis  is 
only  the  beginning  of  what  will  hopefully  *e  a  continuing  effort  here  at  AFIT.  Although  several 
texts  were  referenced  in  this  effort,  one  document  was  extremely  valuable.  This  important  work 
was  Threshold  Signals  by  James  L.  Lawson  and  George  E.  Uhlenbeck. 

I  am  deeply  indebted  to  my  thesis  advisor,  Dr.  Vittal  P.  Pyati,  for  his  help  in  preparing  this 
thesis.  He  was  able  to  keep  me  on  a  consistent  track  when  I  would  begin  to  wander.  I  also  wish  to 
thank  Lt  Col  Meer  and  Lt  Col  Norman  who  both  provided  valuable  assistance  toward  this  effort. 
They  also  sat  on  my  thesis  committee,  for  which  I  am  indebted.  I  have  to  especially  thank  my  wife, 
Sue,  for  allowing  me  the  flexibility  to  work  at  all  hours  at  home  without  interruption.  It  really 
made  a  difference. 


Paul  E.  Bishop 
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.  Abstract 

A  comprehensive  study  of  the  effects  of  noise  jamming  on  various  pulsed  radar  receiver  types  is 
proposed.  First,  a  literature  review  on  the  sources  of  electrical  noise  and  various  jamming  waveforms 
is  presented.  The  power  spectra  for  three  noise  jamming  waveforms  are  rigorously  derived.  These 
cases  are  shot  noise  in  a  parallel-plate  diode,  a  series  of  pulses  with  random  amplitudes  and  intervals, 
and  for  a  series  of  pulses  with  ra..^om  amplitudes,  spacing  and  phases.  Next,  the  probability  density 
functions  for  the  output  of  linear  and  quadratic  detectors  are  developed.  For  each  detector  type,  the 
probability  density  functions  for  the  cases  of  noise  only  and  signal  plus  noise  are  derived.  Finally, 
the  effects  of  CW  interference  on  conventional  pulsed  radar  receivers  is  accomplished.  This  analysis 
shows  that  the  effect  of  a  pure  CW  tone  as  an  interfering  signal  is  that  it  increases  the  minimum 
power  required  for  detection.  As  the  CW  tone  is  detuned  from  the  signal  carrier  frequency  the 
minimum  power  required  for  detection  increases  until  the  detuning  is  so  great  that  the  effects  of 
the  CW  is  negligible.  This  thesis  lays  the  ,  round  work  for  a  much  broader  future  study  of  the 


A  STUDY  OF  THE  EFFECTS  OF  NOISE  JAMMING 


ON  RADAR  RECEIVERS 

L  Introduction 

A  study  of  the  effects  of  noise  jamming  on  radar  receivers  is  proposed.  The  result  of  this  study 
will  be  a  comprehensive  reference  on  the  effects  of  various  types  of  noise  jamming  on  different  types 
of  radar  receivers. 

LI  Background 

Ll.l  Radar  Principles  Radar  is  an  acronym  for  radio  detection  and  ranging.  A  radar 
system  can  measure  the  range  to  a  target  by  transmitting  a  pulse  of  electromagnetic  energy,  and 
measuring  the  time  the  pulse  takes  to  reflect  off  a  target  and  return  to  the  radar  receiver.  The  time 
for  a  pulse  to  travel  one  nautical  mile  and  return  is  called  a  radar  mile.  A  radar  mile  is  defined  as 
12,35 nsec  (6:8).  However,  before  any  information  can  be  obtained  from  a  transmitted  radar  signal, 
that  signal  must  first  be  detected. 

A  radar  receiver  detects  a  target  by  comparing  the  voltage  of  the  received  signal  with  a  preset 
threshold  voltage.  If  the  received  voltage  is  greater  than  the  threshold,  the  radar  is  said  to  have 
detected  a  target.  Radar  receivers  are  designed  to  meet  two  main  opposing  specifications.  These 
are  the  probability  of  detection  (Pd)  and  probability  of  false  alarm  (P/a).  A  radar’s  probability  of 
detection  is  a  measure  of  how  often  the  radar  will  actually  detect  a  target  when  a  target  is  present. 
The  probability  of  false  alarm  is  a  measure  of  how  often  a  radar  will  decide  that  a  target  is  present 
when  there  is  not  a  target  present.  The  parameters  of  Pd  and  P/a  are  used  to  determine  a  radar’s 
required  signal-to-noise  ratio  (SNR)  and  the  minimum  required  signal  power  (18.28).  An  increase 
in  the  noise  power  will  result  in  a  decrease  in  the  signal-to-noise  ratio.  The  decrease  in  SNR,  which 
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is  due  to  an  increase  in  noise  power,  causes  a  decrease  in  the  radar’s  probability  of  detection  and 
an  increase  in  the  radar’s  probability  of  false  alarm. 

i 

LL2  Radar  Receiver  Types  There  are  several  types  of  radar  receivers  in  use  today.  Four  of 
the  most  common  radar  receiver  designs  can  be  grouped  into  four  main  categories.  These  are  the 
conventional  receiver,  the  matched  filter  receiver,  the  pulse  compression  receiver,  and  the  constant 
false  alarm  rate  (CFAR)  receiver. 

The  conventional  receiver  uses  a  a  bilinear  element  followed  by  a  threshold  comparison  device 
to  detect  the  presence  of  a  radar  pulse.  One  example  of  a  conventional  re  eiver  is  a  simple  envelope 
detector,  also  known  as  a  linear  detector.  A  se-  vnd  type  of  conventional  receiver  employs  a  quadratic 
detector.  This  type  of  detector  Is  also  referred  to  as  a  square  law  detector.  A  matched  filter  receiver 
is  designed  to  maximize  ihe  output  signal-to-noise  ratio  (18:369).  The  matched  filter  compares  the 
received  waveform  with  a  stored  replica  of  the  transmuted  pulse.  One  commonly  used  method  of 
implementing  the  matched  filter  is  with  a  correlator  (18:375).  A  correlator  multiplies  a  received 
signal  with  a  reference  signal,  which  is  matched  to  the  transmitted  pulse,  and  then  integrates  the 
resulting  product  over  a  specified  time  period.  The  typ,  of  pulse  compression  recover  analyzed  is 
the  linear  FM  chirp  receiver.  This  type  of  receiver  is  employed  in  high  range  resolution  imaging 
radars.  The  CFAR  receiver  will  be  the  last  receiver  structure  analyzed.  The  CFAR  receiver  employs 
an  adaptive  threshold  to  maintain  a  constant  probability  of  false  ala:  m  (18:39).  These  four  receiver 
types  are  the  most  common  structures  found  in  search  radar  systems. 

1.1.8  Noise  Principles  Electrical  noise  is  electromagnetic  energy  across  a  broad  range  of 
frequencies,  which  is  characterized  by  random  fluctuations  of  amplitude,  frequency,  phase,  etc. 
Noise  is  the  chief  parameter  that  limits  a  radar  receiver’s  sensitivity  (18:18).  Noise  is  also  the 
primary  factor  affecting  a  receiver’s  Pd  and  Pja.  There  are  several  sources  of  electrical  noise.  Two 
main  categories  of  noise  are  the  noise  that  originates  within  the  receiver  and  noise  external  to  the 
receiver.  Thermal  noise  is  the  primary  type  of  noise  generated  in  a  radar  receiver.  The  power  of 
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thermal  noise  is  directly  proportional  to  the  temperature  of  the  receiver.  Thermal  noise  is  present  to 
a  certain  degree  in  all  electrical  devices.  Some  types  of  external  noise  sources  are  man  made  noise, 
cosmic  noise  and  electronic  countermeasures  (ECM).  The  effects  of  thermal  noise  can  be  reduced 
by  proper  leceiver  design,  but  can  never  M  totally  eliminated.  Radar  receivers  are  designed  for  a 
specific  Pd  and  Pja  given  the  known  value  of  thermal  noise  generated  in  the  receiver. 

A  radar  signal  must  travel  through  an  environment  where  electrical  noise  is  present.  A  radar 
receiver  receives  a  signal  that  is  the  summation  of  the  desired  radar  waveform  and  the  external 
noise  (15:36).  A  large  value  of  noise  power  may  cause  the  receiver  to  announce  the  presence  of 
a  target  when  in  reality,  no  target  is  present.  This  false  detection  is  termed  a  false  alarm,  and 
increases  the  radar’s  probability  of  false  alarm.  A  second  effect  of  the  external  noise  is  to  force 
the  receiver  to  conclude  there  is  no  target  when  an  actual  target  is  present,  thereby  decreasing  the 
radar’s  probability  of  detection. 

1.J.4  Noise  Jamming  Noise  jamming  is  one  of  several  active  electronic  countermeasures 
techniques.  The  primary  noise  jamming  method  is  to  produce  a  signal  that  is  as  random,  or  r.oise- 
like,  as  possible  (6:83).  The  optimum  noise  jamming  signal  produced  oy  a  jammer  should  be  a 
close  approximation  of  the  thermal  noise  that  is  produced  by  the  radar  receiver  (20:293). 

One  must  consider  several  variables  when  analyzing  the  noise  jamming  problem.  The  jammer 
power  is  one  of  the  most  critical  parameters.  To  be  effective,  the  level  of  the  jamming  must  be 
greater  than  the  energy  of  the  received  pulsed  radar  signal  (6:89).  The  statistical  characteristics 
of  the  noise  jamming  waveform  are  also  important.  For  a  noise  jamming  signal  to  be  effective,  the 
statistics  of  the  noise  should  be  as  close  to  the  statistics  of  the  receiver  thermal  noise  as  possible. 
Similarly,  a  deception  jammer  attempts  to  mimic  a  radar  target  return.  The  primary  goal  of  noise 
jamming  is  to  increase  the  radar’s  Pja  and  decrease  the  radar’s  Pd  by  producing  a  high  power  noise 
waveform  that  approximates  the  thermal  noise  already  present  in  the  receiver. 
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1.2  Problem 


Adequate  employment  of  radar  systems,  in  an  environment  where  noise  jamming  is  used  by 
an  enemy,  requires  knowledge  of  the  effects  of  the  noise  jamming  on  the  radar  receiver.  Conversely, 
effective  noise  jamming  against  an  enemy  radar  requires  the  knowledge  of  the  effects  of  the  noise 
jamming  on  the  enemy’s  radar  receiver.  Presently,  no  comprehensive  study  exists  on  the  effects  of 
various  noise  jamming  waveforms  on  different  types  of  radar  receivers’  probability  of  false  alarm 
and  probability  of  detection.  A  comprehensive  study  includes  the  characterization  of  the  sources 
of  electrical  noise,  the  statistical  parameters  of  different  noise  jamming  waveforms,  and  an  analysis 
of  the  effects  of  these  jamming  waveforms  on  various  radar  receiver  types.  The  goal  of  this  thesis 
effort  is  to  provide  such  a  comprehensive  study.  The  result  of  this  effort  will  be  a  single  source  of 
reference  for  the  effects  of  noise  jamming  on  radar  receivers. 

1.3  Assumptions 

Several  assumptions  must  be  made  in  the  analysis  of  the  effects  of  noise  jamming  on  radar 
receivers.  The  effects  of  jammer  antenna  gain,  target  characteristics,  range,  and  jammer  power  can 
be  disregarded.  Instead,  the  jammer- to-signal  (J /S)  ratio  will  be  used  in  place  of  these  parameters 
(1:14-11-14-12).  By  the  same  reasoning,  propagation  effects  will  be  disregarded.  The  radar  receiver 
will  be  assumed  to  have  linear  response  characteristics  up  to  the  detector  stage.  This  assumption  is 
a  reasonable  approximation,  which  allows  for  much  simpler  calculations.  However,  the  detector  is 
generally  a  nonlinear  network  and  cannot  be  modeled  with  linear  characteristics.  For  all  analyses, 
only  a  single  pulse  will  be  considered.  For  the  multiple  pulse  case  there  are  simple  techniques 
available  to  determine  the  required  signal-to-noise  ratfo.  A  lossless  system  will  be  assumed.  Any 
system  losses  will  be  a  constant  factor  that  will  not  impact  the  effects  of  the  noise  jamming. 
Finally,  the  radar  receiver  will  be  assumed  not  to ,  incorporate  any  specific  electronic  counter¬ 
countermeasures  (ECCM).  All  these  assumptions  lead  to  a  lossless  receiver  model  that  is  linear  up 
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to  the  detector  stage,  which  is  affected  only  by  the  receiver’s  thermal  noise  and  the  external  noise 
jamming.  Therefore,  the  impact  on  the  modeled  receiver’s  Pd  and  Pja  will  be  due  solely  to  the 
radar  jamming. 

1.4  Approach 

The  study  of  the  effects  of  noise  jamming  will  be  approached  in  thiee  parts.  First,  the 
sources  of  electrical  noise  will  be  characterized.  Next,  the  noise  jamming  waveform  statistics  will 
be  modeled.  Finally,  the  noise  jamming  waveform  statistics  will  be  applied  to  the  four  different 
receiver  models  to  determine  the  effects  of  the  jamming  on  the  receiver’s  Pja  and  Pd . 

The  sources  of  electrical  noise  will  be  characterized  by  their  statistical  properties,  specifically 
by  their  probability  density  and  power  spectral  density.  The  two  primary  sources  of  noise  to  be 
modeled  are  thermal  noise  and  shot  noise.  Both  of  these  sources  can  be  modeled  as  white  (uniform 
power  spectral  density)  and  Gaussian  (normal)  density.  This  model  of  noise  most  closely  resembles 
the  noise  present  in  the  receiver.  Shot  noise  and  thermal  noise  are  the  primary  sources  of  noise 
used  in  radar  noise  jammers.  These  two  noise  sources  are  the  most  popular  because  they  are  simple 
to  generate  by  radar  jammers. 

After  the  sources  of  noise  have  been  modeled,  the  noise  jamming  waveforms  must  be  charac¬ 
terized.  As  with  the  noise  sources,  the  noise  jamming  waveforms  are  modeled  by  their  statistical 
properties.  Although  the  noise  jamming  waveforms  under  consideration  approach  the  ideal  receiver 
thermal  noise  characteristics,  there  are  differences  that  must  be  considered  in  the  analysis  of  their 
effects  on  the  radar  receivers.  These  differences  are  reflected  in  the  noise  power  spectral  density 
and  probability  density  function.  The  noise  jamming  waveforms  to  be  modeled  are  direct  noise 
amplified  (DINA),  frequency  modulated  (FM)  by-noise,  and  amplitude  modulated  (AM)  by-noise. 
Also,  continuous  wave  (single  frequency)  and  random  pulse  waveforms  must  be  considered.  These 
waveforms  are  used  extensively  in  current  radar  noise  jammers. 
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The  analysis  of  the  jamming  effects  on  the  receivers  will  follow  techniques  demonstrated  by 
La»vson  and  Uhlenbeck  (11)  and  those  outlined  by  Bennighof  and  others,  as  well  as  several  other 
radar  text  books  (l).  The  statistical  characteristics  of  the  noise  jamming  waveform  will  be  applied 
to  the  transfer  characteristic  of  the  radar  receiver  to  determine  the  output  noise  characteristics. 
The  result  of  this  analysis  will  be  the  probability  density  and  power  spectral  density  of  the  noise 
output  of  the  receiver.  This  will  be  accomplished  for  the  cases  where  an  actual  target  is  present 
and  when  a  target  is  not  present.  One  technique  to  analyze  the  effects  of  the  radar  jamming  is  to 
apply  the  detectability  criteria  outlined  by  Lawson  and  Uhlenbeck  (11:161-165).  The  application  of 
these  criteria  will  give  an  indication  to  the  required  power  levels  for  detection.  A  second  technique 
is  to  apply  the  two  probability  densities  to  the  known  decision  technique  for  each  radar  receiver 
type.  This  technique  is  similar  to  the  likelihood  ratio  test.  The  results  of  the  likelihood  ratio  will 
determine  the  regions  where  the  Pd  and  Pja  can  be  calculated.  The  resulting  Pd  and  Pja  will  be 
compared  with  the  Pd  and  Pja  for  the  receiver  when  no  noise  jamming  is  present.  The  difference 
in  the  receiver's  Pd  and  Pja  will  be  due  to  the  noise  jamr  .ing.  This  analysis  will  be  accomplished 
for  each  of  the  noise  jamming  waveforms  through  each  of  the  separate  receiver  types. 

1.5  Scope 

This  study  is  limited  to  the  analysis  of  the  effects  of  noise  jamming  on  search  radars.  However, 
the  techniques  used  ar^  applicable  to  the  analysis  of  the  effects  on  tracking  radars  as  well.  Four 
types  of  noise  jamming  waveforms  will  be  used  in  the  analysis.  These  waveforms  are  continuous 
wave  (CYv)  interference,  DINA,  AM  by-noise  and  FM  by-noise.  These  waveforms  are  the  major 
jamming  signals  employed  in  current  noise  jammers  (1.14-9).  Only  pulsed  radars  will  be  considered. 
With  some  slight  differences,  the  methods  employed  in  the  analysis  can  be  aj  plied  to  other  radar 
types,  such  as  continuous  wave  (CW)  and  pulse  doppler  (PD)  radars.  Four  basic  pulsed  radar 
receiver  types  will  be  analyzed  in  this  study.  These  receiver  typos  are  the  conventional  receiver, 
the  matched  filter  receiver,  the  pulse  compression  receiver  and  the  CFAR  receiver. 
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This  thesis  effort,  as  a  part  of  a  much  larger  study,  is  limited  to  the  analysis  of  the  effects  of 
continuous  wave  interference  on  two  types  of  conventional  receiver  detector  structures:  the  linear 
detector,  and  the  quadratic  detector.  The  analysis  is  limited  to  the  application  of  the  deflection 
criterion  outlined  by  Lawson  and  Uhlenbeck  to  determine  the  minimum  signal  power  required  for 
detection  (11:161).  Although  this  effort  is  limited  to  the  analysis  of  CW  interference  on  conventional 
receivers,  the  techniques  developed  here  can  be  easily  applied  to  the  analysis  of  the  effects  on  other 
receiver  structures  by  various  jamming  waveforms. 

L6  Summary 

In  this  chapter,  a  study  of  the  effects  of  noise  jamming  on  various  radar  receivers  was  proposed. 
This  chapter  also  provided  the  necessary  background,  assumptions,  and  the  approach  to  accomplish 
this  study.  This  chapter  also  outlined  the  scope  of  this  particular  thesis  effort.  Chapter  II  is  a 
literature  search  on  the  sources  of  electrical  noise  and  on  various  noise  jamming  waveforms.  In 
Chapter  III  the  power  spectra  for  three  specific  jamming  waveforms  are  derived.  The  probability 
density  functions  for  the  output  of  linear  and  quadratic  detectors  are  developed  in  Chapter  IV. 
Probability  density  functions  for  the  cases  of  noise  only  and  of  signal  plus  noise  are  derived  for  each 
detector  type.  The  effects  of  continuous  wave  interference  on  linear  and  quadratic  detectors  are 
analyzed  in  Chapter  V.  Lastly,  Chapter  VI  contains  conclusions  and  recommendations  for  further 
work  in  this  study. 
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II.  Literature  Review 


Before  analyzing  the  effects  of  noise  jamming  on  radar  receivers,  one  needs  an  understanding 
of  two  main  subject  areas.  These  subject  areas  are  the  sources  of  electrical  noise  used  in  radar 
jammers  and  noise  jamming  waveforms.  The  critical  parameters  to  consider  for  the  sources  of  noise 
are  the  signal’s  probability  density  function  and  power  spectral  density.  These  two  parameters 
are  also  important  to  the  study  of  noise  jamming  waveforms.  A  great  deal  of  research  has  be^n 
conducted  in  these  arsas.  Therefore,  the  goal  of  this  literature  review  is  to  determine  the  most 
recent  research  in  these  areas,  and  the  results  of  that  research. 

2.1  Noise  Sources 

The  basis  of  most  radar  jamming  signals  is  a  source  of  electrical  noise.  Any  noise  source 
may  be  generally  categorized  in  one  of  two  ways,  either  by  its  power  spectral  density  (PSD)  or 
the  probability  density  function  (PDF)  of  one  of  its  parameters  (usually  the  amplitude  of  the  noise 
voltage).  The  PDF  of  a  signal  is  a  measure  of  how  likely  a  particular  amplitude  of  the  signal  is, 
whereas  the  PSD  of  a  signal  is  a  measure  of  how  much  power  is  contained  in  specific  frequency 
bands.  Turner  and  others  point  out  that  “  . . .  the  ECM  goal  is  to  produce  white  Gaussian  noise  in 
the  victim  radar  receiver”  (19:118).  Gaussian  refers  to  the  PDF  of  the  amplitude  of  the  noise.  This 
type  of  noise  has  a  Gaussian  or  normal  amplitude  PDF  (represented  by  the  familiar  bell  shaped 
curve).  White  noise  is  noise  that  has  a  uniform  or  flat  PSD  over  all  frequencies,  analogous  to  the 
concept  of  white  light,  which  contains  light  at  all  visible  wavelengths  (9.48).  The  true  goal  of  noise 
jamming  is  to  produce  a  Rayleigh  distributed  noise  output  from  the  victim  receiver’s  detector(7.66). 
The  Rayleigh  PDF  is  closely  related  to  the  Gaussian  PDF.  Rayleigh  distributed  noise  describes  the 
envelope  of  the  output  of  a  radar’s  intermediate  frequency  (IF)  filter  when  the  input  of  the  filter  is 
Gaussian  noise  (18:24).  It  is  sufficient  for  the  analysis  of  the  effects  of  noise  on  radar  receivers  that 
the  PSD  of  the  noise  be  uniform  only  over  a  specific  frequency  bandwidth  of  interest,  generally  the 
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receiver  bandwidth.  There  are  several  sources  of  white  Gaussian  noise.  The  most  common  sources 


are  thermal  noise  and  shot  noise. 

Thermal  noise  is  due  to  the  random  motion  or  fluctuations  of  electrons  in  conductors  (9:47). 
When  quantum  mechanical  effects  are  neglected,  this  type  of  noise  has  a  flat  or  uniform  power 
spectral  density,  directly  proportional  to  the  temperature  of  the  device  (21:9).  By  the  central  limit 
theorem,  which  states  that  for  the  sum  of  a  large  number  of  independent  samples  the  limiting 
form  of  the  PDF  is  Gaussian,  thermal  noise  has  a  Gaussian  PDF  with  a  mean  value  of  zero  and  a 
variance  equal  to  the  mean  square  value  of  the  noise  voltage  (9:48-49).  Thermal  noise  is  present  to 
a  certain  extent  in  all  electrical  devices  and  can  be  exploited  as  a  noise  source  for  radar  jammers. 
Another  source  of  electrical  noise  commonly  used  in  radar  jammers  is  shot  noise. 

Shot  noise  results  from  the  flow  of  electrons  in  active  devices  such  as  vacuum  tubes  and 
semiconductors.  In  vacuum  tubes,  the  shot  noise  is  caused  by  the  current  induced  from  the  ran¬ 
dom  emission  of  electrons  from  the  heated  cathode  (4:112).  The  shot  noise  in  semiconductors  is 
caused  by  currei  t  flows  generated  when  carriers  (electrons  and  holes)  randomly  crossing  the  p-n 
barrier  (22:93).  Davenport  and  Root  provide  the  development  of  the  PSD  and  PDF  of  shot  noise 
in  thermionic  vacuum  tubes  for  both  the  temperature-limited  diode  and  the  space-charge-limited 
diode  cases  (4:112-143).  They  showed  that  for  frequencies  less  than  the  reciprocal  of  the  electron 
transit  time  in  the  diode,  the  PSD  of  shot  noise  can  be  approximated  by  a  flat  spectrum  (4:123). 
This  spectrum  is  proportional  to  the  average  current  through  the  diode,  as  opposed  to  the  device 
temperature  as  is  the  case  for  thermal  noise.  Shot  noise  is  also  a  phenomenon  present  in  semicon¬ 
ductors.  Van  der  Ziel  developed  expressions  for  shot  noise  in  semiconductor  diodes  for  both  the 
low  frequency  and  high  frequency  cases  (22:93-109).  He  also  presents  an  expression  for  the  PSD  of 
noise  from  a  back-biased  p-n  diode,  and  states  that  large  amounts  of  noise  power  can  be  generated 
in  this  case  (22:48).  Penney  and  others  describe  the  back-  biased  diode  as  an  “extremely  hot  noise 
generator”  (9:71).  By  the  application  of  the  central  limit  theorem,  the  PDF  of  shot  noise  in  both 
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vacuum  tube  diodes  and  semiconductor  diodes  can  be  shown  to  approach  a  Gaussian  density.  Shot 
noise,  like  thermal  noise,  is  one  excellent  source  of  white  Gaussian  noise. 

The  properties  of  Gaussian  noise  make  it  an  ideal  choice  for  radar  jamming  signals.  Turner 
and  others  define  a  measure  of  jamming  performance,  the  noise  quality  (19:117-122).  The  noise 
quality  is  a  measure  of  how  closely  the  PDF  of  the  jamming  signal  approaches  that  of  a  Gaussian 
density.  The  authors  show  that  a  high  value  of  noise  quality  can  lead  to  a  reduction  in  the  required 
jamming-to-signal  ratio  (J/S)  (19:117).  For  a  given  level  of  effectiveness  against  a  pulsed  radar, 
this  reduction  in  J/S  can  be  as  much  as  17  dB  (a  factor  of  50)  over  a  low  noise  quality  signal 
that  produces  the  same  effect  (19:117).  Based  on  the  work  done  by  Turner  and  others,  Knorr 
and  Karantanas  outline  a  computer  simulation,  in  use  at  the  Naval  Postgraduate  School,  for  the 
optimization  of  various  jamming  waveforms  (10:273-277).  This  computer  simulation  compares  the 
PDF  of  the  jamming  waveform  with  a  Gaussian  PDF,  having  the  same  mean  and  variance,  to 
arrive  at  a  value  for  noise  quality,  as  defined  by  Turner  and  others  (10:274).  Work  continues  at  the 
Naval  Postgraduate  School  in  this  area.  Even  though  Gaussian  noise  is  the  most  common  source  of 
noise  for  radar  jammers,  other  noise  densities  can  also  be  used.  Two  such  noise  sources  are  random 
pulses  and  single  frequency,  or  continuous  wave  (CW)  noise. 

Random  pulses  can  be  used  as  one  source  of  noise  for  a  jamming  signal.  Lawson  and  Uhlenbeck 
present  a  thorough  analysis  of  the  power  spectra  of  random  pulse  signals  (11:43-46).  They  provide 
the  derivation  of  the  power  spectra  for  a  series  of  pulses  with  random  amplitudes,  random  pulse 
repetition  intervals,  and  random  phases  (11:43-46).  The  power  spectrum  for  a  series  of  pulses 
with  random  amplitudes  is  a  continuous  spectrum  that  has  the  same  shape  as  the  power  spectrum 
of  a  single  pulse  (11:44).  However,  the  power  spectrum  for  a  series  of  pulses  with  random  pulse 
repetition  intervals  is  dependent  on  both  the  spectrum  of  a  single  pulse  and  the  probability  density 
of  the  pulse  intervals  (11:44).  Maksimov  and  others  present  expressions  for  the  average  values  of 
pulse  duration  and  pulse  spacing  for  random  pulse  waveforms  (12:44).  Random  pulse  waveforms 
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may  be  used  alone  as  a  jamming  source.  A  series  of  random  pulses  can  also  be  used  to  modulate 
other  noise  waveforms. 

A  continuous  wave  signal  (single  frequency  tone)  can  also  be  a  source  of  interference  to  a  radar 
receiver.  For  any  given  frequency,  the  CW  tone  can  have  a  randomly  varying  amplitude  and  phase. 
Jordan  and  Penney  present  an  expression  for  the  PDF  of  the  voltage  of  the  CW  interference  (9:46- 
47).  The  PDF  of  the  amplitude  of  the  CW  interference  iollows  a  somewhat  parabolic  distribution 
between  the  peak  voltage  values  (9:46-47).  For  the  density  of  the  phase  of  CW  interference,  in  the 
case  of  a  pulsed  radar  it  is  equally  likely  that  the  CW  interference  is  in  phase  with  the  received 
radar  pulse  as  it  is  out  of  phase  with  the  pulse  (11.336).  As  is  the  case  for  random  pulses,  CW 
tones  can  be  used  alone  as  a  jamming  signal. 

2.2  Jamming  Waveforms 

There  are  several  types  of  waveforms  employed  in  noise  jammers.  This  section  will  discuss 
three  main  types  of  barrage  jamming  waveforms.  Barrage  jamming  refers  to  radar  jamming  over  a 
broad  band  of  frequencies,  as  opposed  to  spot  jamming  which  is  defined  as  the  jamming  of  a  specific 
frequency  (6:253,268).  These  barrage  jamming  waveforms  are  direct  noise  amplified  (DINA),  fre¬ 
quency  modulated  (FM)  by-noise  and  amplitude  modulated  (AM)  by-noise.  The  majority  of  what 
follows  in  this  section  is  taken  from  a  chapter  of  Electronic  Countermeasures  entitled  “Effectiveness 
of  Jamming  Signals”  by  Bennighof  and  others  (1:14-1-14-65). 

DINA  is  merely  amplified,  bandlimited  Gaussian  noise  (1:14-9).  Over  the  passband  of  interest, 
DINA  can  be  considered  white  (uniform  PSD)  (1:14-9).  DINA  can  be  easily  generated  by  amplifying 
low  level  thermal  noise  (1:14-9-14-10).  DINA  has  all  of  the  properties  of  white  Gaussian  noise. 
For  the  purpose  of  analyzing  the  effects  of  DINA  on  radar  receivers,  the  noise  can  be  considered 
to  be  additive.  That  is,  the  noise  is  added  to  the  radar  signal  prior  to  reception.  The  analysis  of 
the  effects  of  additive  white  Gaussian  noise  (AWGN)  is  presented  in  several  radar  texts.  Skolnik 
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gives  an  analysis  of  the  effects  of  AWGN  on  a  radar’s  probability  of  detection  and  probability  of 
false  alarm  (18:23-29).  There  are  other  jamming  waveforms  employed  in  radar  jammers  that  are 
not  truly  Gaussian,  but  approach  a  true  white  Gaussian  noise  source. 

In  addition  to  DINA,  FM-by-noise  is  one  frequently  employed  class  of  noise  jamming  wave¬ 
forms.  There  are  two  distinct  types  of  FM-by-noise:  FM-by-WB  (wideband)  noise,  and  FM-by-LF 
(low-frequency)  noise.  (1:14-10).  Each  type  of  FM-by-noise  waveform  has  different  effects  on  radar 
receivers.  FM-by-noise  is  used  extensively  in  current  radar  jammers. 

One  frequently  used  noise  jamming  waveform  is  FM-by-WB  noise.  FM-by-WB  noise  produces 
a  spectrum  nearly  the  same  as  that  of  DINA.  (1:14-10).  Bennighof  and  others  showed  that  the  noise 
output  of  a  radar  receiver’s  IF  filter  due  to  FM-by-WB  noise  can  be  the  same  as  the  IF  filter  output 
due  to  DINA  (1:14-10).  The  effects  of  FM-by-WB  noise  are  the  same  as  those  of  DINA  when  the 
bandwidth  of  the  barrage  jamming  waveform  is  greater  than  the  noise  bandwidth,  and  greater  than 
the  receiver  bandwidth  (1:14-22).  Another  FM-by-noise  waveform  used  in  several  noise  jammers  is 
FM-by-LF  noise. 

FM-by-LF  noise  is  produced  in  much  the  same  way  as  FM-by-WB  noise,  but  the  noise  signal  is 
confined  to  a  much  narrower  modulating  bandwidth,  a  bandwidth  less  than  the  victim  radar’s  band¬ 
width  (1:14-10).  Bennighof  and  others  point  out  that  FM-by-LF  noise  has  two  main  advantages: 
FM-by-LF  noise  produces  more  power  out  of  the  victim  radar’s  IF  filter  than  does  FM-by-WB 
noise,  and  FM-by-LF  noise  is  more  effective  in  the  jamming  of  radar  systems  which  use  plan  po¬ 
sition  indicator  (PPI)  displays  (1:14-10).  They  also  point  out  that  FM-by-LF  noise  has  one  main 
disadvantage,  that  its  effects  are  easily  negated  by  electronic  counter-countermeasures  (ECCM) 
(1:14-10). 

FM-by-noise  is  an  extremely  popular  jamming  waveform.  Cassara  and  others  describe  a 
technique  to  generate  a  uniform  PSD  jamming  signal  from  a  Gaussian  source  using  an  FM-by-noise 
method.  (2:330).  The  technique  they  present  is  to  pass  Gaussian  noise  through  a  nonlinear  network, 
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followed  by  a  fr.q:ftncy  modulator  (2:330-332).  The  result  of  their  work  was  an  extremely  good 
source  of  uniform  PSD  noise  over  an  electronically  controlled  bandwidth  (2:332).  The  analysis 
presented  by  Knorr  and  Karantanas  was  mainly  geared  towards  the  study  and  optimization  of 
FM-by-noise  waveforms  (10:273-277).  The  results  of  their  study  gave  an  indication  of  methods  to 
improve  the  noise  quality  of  an  FM-by-noise  waveform  (10:275-277). 

Another  jamming  waveform,  although  used  much  less  frequently  than  FM-by-noise,  is  AM-by- 
noise.  The  AM-by-noise  effects  are  similar  to  those  of  DINA,  but  it  is  difficult  to  produce  wideband 
noise  in  this  manner  (1:14-11).  This  difficulty  arises  because  the  bandwidth  of  the  AM-by-noise 
waveform  is  limited  to  twice  the  bandwidth  of  the  modulating  noise  (1:14-11).  This  bandwidth 
limitation  results  from  the  property  that  the  bandwidth  of  any  AM  (double  sideband)  signal  is 
equal  to  twice  the  bandwidth  of  the  modulating  signal.  Another  reason  that  AM-by-noise  is  not  as 
popular  for  use  in  jammers  is  that  part  of  the  total  power  transmitted  by  an  AM-by-noise  jammer 
is  lost  in  the  carrier  (assuming  large-carrier  AM).  In  fact,  no  more  than  50  percent  of  the  total 
power  will  be  in  the  jamming  portion  of  the  signal.  The  remaining  power  is  in  the  carrier  (12:37). 

A  further  jamming  technique  is  one  that  employs  random  pulse  modulation.  With  this 
method,  either  an  FM-by-WB  noise  or  DINA  waveform  is  pulse  modulated.  (1:14-11).  The  re¬ 
sulting  random  pulse  jamming  waveform  employs  random  duration,  spacing,  and  amplitude  of  the 
pulses,  where  the  average  pulse  duration  should  be  the  same  as  the  victim  radar’s  pulse  width 
(1:14-11).  The  effects  of  the  random  pulse  jamming  waveform  are  similar  to  those  produced  by 
the  FM-by-LF  noise  waveform  (1:14-11).  The  most  effective  random  pulse  jamming  signal  is  one 
in  which  the  amplitude,  pulse  duration  and  pulse  interval  are  all  varied  randomly.  However,  such 
a  signal  is  difficult  to  generate  (12:42).  It  is  much  simpler  in  practice  to  generate  a  signal  with 
constant  amplitude  pulses  where  the  duration  and  pulse  interval  are  varied  randomly  (12:42). 
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2.S  Summary 


As  this  chapter  has  pointed  out,  a  great  deal  of  work  has  been  accomplished  in  the  area  of 
noise  jamming.  This  literature  search  has  concentrated  on  the  sources  of  electrical  noise  and  on 
noise  jamming  waveforms.  Equations  for  the  PDF’s  and  PFD’s  of  several  different  noise  sources 
have  been  developed.  There  has  also  been  much  effort  in  the  area  of  analyzing  noise  jamming 
waveforms.  The  goal  of  this  review  was  to  determine  the  most  current  information  on  these  two 
subjects. 
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III.  Power  Spectra 


One  of  the  main  parameters  to  consider  when  studying  the  effects  of  jamming  on  radar 
receivers  is  the  power  spectral  density  (PSD)  of  the  noise  waveform.  The  PSD  describes  the 
frequency  distribution  of  power  in  the  noise  jamming  waveform.  In  this  chapter,  the  power  spectra 
for  three  cases  is  developed.  First,  the  PSD  for  the  shot  noise  in  a  parallel-plate  diode  is  derived. 
Next,  the  PSD  for  a  series  of  pulses  of  random  amplitudes  and  spacing  is  developed.  From  this  PSD, 
the  special  cases  of  the  PSD  for  random  binary  transmission  and  the  boxcar  spectrum  are  derived. 
Lastly,  the  power  spectrum  for  the  general  case  of  a  series  of  pulses  with  random  amplitudes, 
spacing,  and  phases  is  derived.  For  this  general  case,  the  special  cases  of  the  PSD  of  the  random 
binary  transmission  and  of  the  random  telegraph  signal  are  developed. 

S.l  Shot  Noise 

In  this  section,  the  power  spectral  density  of  the  shot  noise  in  a  temperature-limited  parallel- 
plate  thermionic  vacuum  tube  diode  is  be  derived.  Much  of  the  development  here  follows  the 
treatment  of  this  subject  by  Davenport  and  Root  (4:112-124). 

The  anode  current  pulse  due  to  the  travel  of  a  single  electron  is  represented  by  the  following 
equation: 


ie(t)  = 


l  *(£) 

1° 


0  <  t  <  ra 
elsewhere 


(3.1) 


where 

q  —  charge  per  electron  (  -1.6  x  10~19  Coulomb) 
ra  =  transit  time 
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The  total  current  flowing  through  the  vacuum  tube  diode  is  the  sum  of  the  individual  electron 
current  pulses.  The  total  current  in  a  time  interval  (— T,  T )  resulting  from  the  flow  of  K  electrons 
is  represented  by  the  following  equation: 


K 


*^(0  —  ^  (*  tk) 
*= 1 


for  - T<t<T 


(3.2) 


where 

tk  =  emission  time  of  the  kih  electron 

The  number  of  electrons  emitted  in  a  given  time  interval  is  assumed  to  have  a  Poisson  dis¬ 
tribution.  The  probability  of  I(  electrons  emitted  in  a  time  interval  r  is  given  by  the  following 
expression  for  the  Poisson  distribution: 

P(K,r)=  (-gT>  (3.3) 

where 

n  =  average  number  of  electrons 

The  following  statistical  properties  are  obtained  from  the  Poisson  distribution: 


X- - -V 

II 

S? 

(3.4) 

Et  (K2)  =  nr  +  (nr)2 

(3.5) 

Var r(K)  =  nr 

(3.6; 
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8.1.1  Average  Current  The  total  diode  current  I(t)  is  a  function  of  tlu  ^al  number  of 
electrons  emitted,  K}  and  the  emissi  a  times  of  each  electron,  the  V s.  These  terms  are  all  random 
variables.  The  average  current  is  obtained  from  the  expected  value  of  the  total  current  I(t): 


E[I(t)}=  r  ...  r  I(t)p ttK>K)dh...dtKdK  (3.7) 

J —oo  J- 00 

The  joint  PDF  of  the  number  of  electrons  and  the  emission  times  may  be  written 


p(h . tK>K)  =  p(tu...,tK  I  K)p{K )  (3.8) 

The  emission  times,  the  Vs,  are  independent,  uniformly  distributed  random  variables  with  proba¬ 
bility  density  function  (4:117-119) 


P(h)  =  < 
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t  <  tk  <  i  +  ?T 
elsewhere 
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Since  the  Vs  are  statistically  independent  random  variables,  the  joint  PDF  of  the  emission  times 
and  the  number  of  electrons  becomes 


P  (*i . tx,K)  =  p  {h ,  |  K) . .  .p  (ti( ,  1 1()  p(I<) 


and  the  expected  value  of  the  total  current  can  be  written 


(3.10) 


m*)) 


rCO  r  CO  rOO 

I  p{K)dK  .../  7(0 

J -oo  J  —  oo  J -oo 


dl\  tuft 

2T  ■"  2T 


/co  roo  /'oo  K 

p(I()dK  ...  'Zicit-tk) 

■CO  J- CO  J  —  OO  =  1 


dt\  dtj( 

2T  2T 
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/oo  ft  /  «  \  k  y»T  />cr  /r 

^ p(/o<w  E (2?)  J Tdti"-J T dtK~i j T ie v ~ **) dtK 

=  £^JKf|:(^)'t(2T)K',?l 


(3.11) 


Let  7(2)  =  E  [!(<)].  The  expected  value  of  the  current  (the  DC  current)  is 


7  -  fW  tjf 

=  wCKmdK 

=  |;f?2T(/0 

=  nq  (3.12) 

8.1.2  Power  Spectral  Density  By  the  application  of  the  Wiener-Khinchine  relationship, 
which  states  that  the  power  spectral  density  of  a  stationary  process  is  obtained  from  the  Fourier 
transform  of  the  autocorrelation  function,  the  power  spectral  density  is  obtained  (17:145).  The 
autocorrelation  function  is  derived  as  follows: 


.dtkdK 


Rj{r)  =  E[I(t)I(t  +  r)] 

/OO  rOO 

•••  /  I(t)I(t  +  T)p(tl}...,tK,K)dti...( 

*00  J  —  OO 

/OO  rOO  K  ft 

•••  /  +  l-'-dticdK  (3.13) 

•OO  J- OO  j  =  1 


Because  the  7^’s  are  independent  and  uniformly  distributed  with  the  probability  density  function 
of  Eq  (  3.9)  the  equation  for  the  autocorrelation  function  may  be  written 
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Ri(t)  =  ["  p(K)dK  f  ■■■  jT  + 

*/  — oo  J  — T  J  —T  ^  ^*  =  1 

-  />>- 


2T 


(3.14) 


There  are  if2  terms  under  the  double  summation:  K  terms  when  i  =  j>  and  K2  -  K  terms 
when  i  ^  j. 

For  i  =  j 


For  i  ^  j 


dt\ 
2 T 


•  /  *e(*  —  <*)*«(<  +  T  —  <i) 

J-T 


dtjf 
2  T 


{&)  Wf-'  J\(t)h(t  +  T)dt 
i  yT 

2T  J  _t‘e(t)*«(<  +  r)d<  (3-15) 


dii 

2T 


•/. 


/p  , 


[±)K  Wf-* 

fT 

x  /  ic{t  +  r  -  tj)dtj 

J-T 

(2?)  S.T>’m 

«’  (±fj  (3.1C) 


Combining  the  above  two  cases,  the  autocorrelation  function  of  the  total  current  is  now 
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p(I<)dK 


-  ci(K2-K)^y 

/  n  \  2  Z*00 

=  Of)  j_J^-K)p(K)iK 

+  [55;  f°  +  T)it 1  J“  K,(K)JK 

/  n  \  2  1  fCO 

=  (^)  ^2T  (/C2  -  K)  +  —E2T  (JO  J  "  ie(t)ie(t  +  r)<tt  (3.17) 


Using  Eqs  (  3.4)  and  (  3.5)  the  expectations  are 


E2t(I<)  =  2Tn 


(3.18) 


E2T  {I<2  -  K)  =  (2Tn)2 


(3.19) 


With  these  expectations,  the  autocorrelation  function  is 


ri(j) 


(^)2  (2^)2  +  -  J  (2TB)  *■,(<)£.(<  +  r)d* 

y*00 

(qn)2  +  n  ic(t)ie(t  +  r)dt 


(3.20) 


The  shot  noise  current  of  interest  is  the  AC  portion  of  the  total  current  (the  fluctuation  about 
the  mean  value).  Let  i(t )  be  the  AC  diode  current: 


(3.21) 


The  autocorrelation  function  of  the  AC  current  is 
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(3.22) 


ie(t)ie(t  +  r)dt 


Since  the  shot  noise  is  a  stationary  process,  the  Fourier  transform  of  the 
function  of  the  process  yields  the  power  spectral  density  of  the  shot  noise: 


StU)  =  ^  Rt(r)e-MT  dr 

/CO  A  CO 

n  /  +  r)dte~*2*JTdT 

•CO  j -CO 


Let  r  =  tf  —  iy  so  the  power  spectral  density  now  becomes 


/CO  f  CO 

ie{t)e-^^'-^dt  /  it(t')dt' 

•CO  J -00 

/CO  fOO 

ic(t)ei2*Jtdt  /  dt' 

■00  */-oo 


Let  G(/)  be  the  Fourier  transform  of  the  current  pulse: 


/oo 

ie(i)e-iwtdt 

■00 


where 


a;  =  2tt/ 


so  that  the  power  spectral  density  is 


S,-(w)  =  n|G(w)|2 


autocorrelation 


(3.23) 


(3.24) 


(3.25) 


(3.26) 
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For  low  frequency  values  (/  4C  l/ra)  the  PSD  is  approximately 


Si(w)  «  n  |G(0)|2  (3.27) 

From  Eqs  (  3.1)  and  (  3.25)  G(0)  =  q  and  it  follows  from  Eq  (  3.12)  that  the  power  spectral  density 
is  approximated  by 


Si(w)  ~  Iq 


(3.28) 


which  is  the  Schottky  formula. 

3J.3  The  Power  Spectral  Density  for  the  Parallel-Plane  Diode  In  this  section  the  power 
spectral  density  for  the  shot  noise  from  the  parallel-plane  diode  will  be  developed.  The  expression 
for  the  power  spectral  density  is  given  by  Eq  (  3.26).  The  current  pulse  is  defined  by  Eq  (  3.1)  and 
the  Fourier  transform  of  the  current  pulse  is 


G(u)  =  %te~iwtdt 
Jo  rl 

=  7^72  "  e_,’WTa  "  3“rae-iWT*  ) 

(iora) 

=  - 2  [0  ”  coso;ra  -  o>ra  sino;ra)  4-  j  (sinwra  -  corQ  coso;ra)]  (3.29) 


The  square  of  the  magnitude  of  the  Fourier  transform  of  the  current  pulse  is 


|G(w)|2  =  —  " j  [(3  ~  cos -  ura  sin  wr,)2  +  (sin  wr0  -  u>ra  coswra)2j 


(wra)4 

{"raf 


|(wra)2  +  2  (1  —  cos  wr„  —  wr„  sin  wr„)  j 


(3.30) 
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Substituting  the  above  equation  into  Eq  (  3.26)  the  power  spectral  density  becomes 


Si( w)  =  --—r  f(u>ra)2  +  2  (1  -  cos UTa  -  wrfl  sinwra)l  (3.31) 

(wra)  »■  J 

From  Eq  (  3.12)  liq  ~  I  and  the  power  spectral  density  becomes 


■S{(w)  =  ~~~r  [(wr0)2  +  2  (1  -  cos UTa  -  wra  sinwr0)] 


(wr0) 


which  is  illustrated  in  Figure  3.1. 


(3.32) 


Figure  3.1.  Parallel-Plane  Diode  Shot  Noise  PSD 


This  equation  agrees  with  the  expressions  for  the  PSD  derived  by  Davenport  and  Root  (4:124)  and 
Papoulis  (14:357-360). 
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8.2  Pulses  with  Random  Amplitudes  and  Intervals 


The  power  spectral  density  for  a  series  of  pulses  with  random  amplitudes  and  pulse  intervals 
is  developed  here.  An  infinite  series  of  pulses  with  random  amplitudes  and  intervals  is  defined  by 
the  following  equation: 


CO 

y(t)=  £  akg(t-kT9-ek)  (3.33) 

k——oo 

where  a *  is  the  pulse  amplitude  with  PDF  p(a),  T0  is  the  average  pulse  repetition  interval,  and  e* 
is  the  deviation  from  T0  of  the  ki\\  pulse  with  PDF  p(e)  and  E(e)  =  0.  The  series  is  truncated  to 
a  finite  interval  from  -NT0  to  +NT0: 


N 

J/w(0  =  13  ak9{i  -  kT0  -  ek)  (3.34) 

k=-N 

A  portion  of  a  series  of  pulses  with  random  amplitudes  and  spacing  is  illustrated  in  Figure  3.2. 
The  Fourier  transform  of  the  sequence  is 


N 

YN(f)  =  G(f)  £  ake-J*'/VT.+tt.)  (3.35) 

where  the  Fourier  transform  of  a  single  pulse  is 


/CO 

g^e-^dt  (3.36) 

•CO 

8.2.1  Power  Spectral  Density  The  expected  value  of  the  power  spectral  density  is  derived 
here.  The  expected  value  will  be  used  since  the  a*  and  terms  are  random  variables.  The  average 
PSD  is  defined  by 
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Figure  3.2.  Series  of  Pulses  with  Random  Amplitudes  and  Spacing 


^«  =  »lm»(S7Ti)5’|y''(/)|2  (3'37) 

where  there  are  approximately  2AT+  1  pulses  between  —NT0  and  4-NTo.  The  square  of  the  mag¬ 
nitude  of  the  Fourier  transform  of  the  series  of  pulses  is  represented  by  the  following  equation: 

\YN(f)\2  =  |G(/)|2  E  akale-j2*WT°+ek'>e+>2*fVT°+ei'>  (3.38) 

k=-Nl=-N 

The  expected  value  of  the  above  equation  is 


!W)I2 


N  N 


|G(/)|2  Xj  J2  [e~*2”fkT°e'i'j2*flT° 

k=-Nl=-N 


x  E  (ckaie-^^e+i2*'1')} 


(3.39) 
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Since  the  ajt’s  and  ej^’s  are  statistically  independent  the  expectation  in  the  above  expression  has 
the  following  values 


E(akate->2**€ke+J2*tc') 


t 

E{cr)  ,  k  =  / 

£2(a)l<M/)|2  . 


(3.40) 


where 


*«(/)  =  p(ey2*J‘dc 

J -CO 


(3.41) 


which  is  the  characteristic  function  of  the  deviation,  e. 

There  are  2AT-f-l  terms  where  k  =  /.  Therefore  the  expression  for  the  square  of  the  magnitude 
of  the  Fourier  transform  becomes 


\YN(f)\2  =  |G(/)|2 


(2 N  +  1  )E(a2)  +  E2(a)  |<h£(/)|2  £  ]£  s-WtT°e+i2*J,T° 

Jt? n 


(3.42) 


The  double  summation  in  the  above  equation  may  be  written 


N  N 

=  —(2N  +  1)  +  V)  V  e-jWkT0f,+j2njlT<> 

kjtl  k=-Nl=-N 


=  -(2N  +  1)  + 


E 

*=-iV 


e>-j2nfkT0 


(3.43) 


so  that  the  expected  value  of  the  magnitude  of  the  square  of  the  Fourier  transform  now  becomes 
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IW)I2 


|G(/)|2  { (2 N  +  1)  [E(a 2)  -  E\a)  |$£(/)|: 


+E2(a)  |$c(/)|2 


AT 


^  e-j2xJkT0 


Jfess-AT 


} 


(3.44) 


The  average  power  spectral  density  is  arrived  at  by  taking  the  limit  of  the  above  expression  as  N 
approaches  infinity: 


Mf)  “  ApSrTi)E|G(/)|! ’{(2W+1)K!)-®2(»)l*<(fl|!l 


+E2(a)  |$.(/)|2 
l 


N 


^  e-j2xJkT, 


k=-N 

2  f  nV  _2\ 


n 


T„ 


\G(f)\UE(a*)-EJ(a)\$c(f)V 


+E2(a )  |$£(/)|2  Jim 


1 


jV— >00  2 N  4*  1 


N 


NP  e-j2*JkT0 
k=-N 


(3.45) 


The  summation  in  the  above  equation  may  be  expanded  using  Euler’s  formula  and  can  be  written 


N 


N 


-1 


e-j2 xJkT0 

kzz-N 


=  i  +  Y^j27rfkTo+  £  e~j2*fkTo 

k= 1  k=-N 

N  N 

_  1  4-  e-j2jr/*T0  _J_  e+j2*JkTo 

k-X  k=l 

N 

—  J  .j.  ^  ^e“i2ff/^ro  ^  e+j2xfkT0^ 

k= 1 
N 

=  1  4-  2  ^  cos  2irfkT0 


(3.46) 
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The  following  expression  for  the  sum  of  cosines  is  given  by  Eq  (.1.342,2)  from  Gradshteyn  and 
Ryzhik  (8:30): 


n 


y  cos  kx 
0 


n  +  1  .  nx  x 

cos  — - — a:  sm  — cosec-  +  1 
Z  z  z 


which  can  be  manipulated  to  the  form 


(3.47) 


N 


y  cos  kx  = 


jV-H  •  Nx 

cos  ~ 


*=i 


(3.48) 


Using  Eq  (  3.48)  the  summation  in  Eq  (  3.46)  becomes 


N 


Y  e-j2*tkT* 

k=-N 


2cos(£^)27rfkT0  sin  N~fkT°- 
1  +  sin  2^~'n 


Using  the  trigonometric  identity 


(3.49) 


cos  A  sin  D  =  ^  sin  ( A  +  B)  —  i  sin  ( A  -  B) 
z  z 


(3.50) 


the  summation  can  be  expressed  as 


V'  fi-j2*fkT0  _  !  ,  sin(2A^  +  l)7r/r0  -sin7r/T0 
~  sin  ~fT0 

_  sin  7r fT0  +  sin  (2 N  +  1  )*fT0  -  sin  tt fT0 
sin  7T  fT0 

_  sin  (2  N_  +  l)vfT0 
sin  ?r  fT0 


(3.51) 


The  limit  now  becomes 
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(3.52) 


1 


lim 

N-+  co  2N  4- 1 


N 

2-  lim  1 

’sin  (2iV  +  l)irfT0' 

JLi  e 
k=-N 

llill  r 

N-*<x>  2N  +  1 

sin  nfT0 

1  2 


By  letting 


ft  =  (2N  +  1) 


a  =  irfTo 


the  limit  can  be  expressed 


lim  — r 
N-+oo  2N  +  1 


N 


e-j2xfkT0 

jl4-at 


=  lim  — 

O— ►oo  ft 


Cl  ( sinafl\2l  ( 7ra2fi\ 

»(— j  <3-53> 


Using  the  following  expression  for  the  limit  of  the  square  of  the  sinc-over-argument  in  the  above 
equation 


lim  2(Wiy 

n-^co  t r  \  aft  /  w 


(3.54) 


leads  to  the  following  form 


1 


lim 

jv-oo  2N  4- 1 


AT 


Y  Wo 


k=-N 


7TCT 

sin2a 


tt5(cv) 


6{a) 


O'* 


sin  * a 


a= o 


(3.55) 


By  L’Hopital’s  rule,  the  following  result  is  achieved: 


cr 


lim  . 

Of-*0  Sill  lct 


=  1 


(3.56) 


which  leads  to  the  following  equation: 
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lim  —7 — - 

N  — ♦  00  2  N  +  1 


N 


^  e-;27r/fcTe 


=  n8(nfT0) 


(3.57) 


Using  the  following  property  of  the  delta  function 


«(«/) = jjj«(/) 


(3.58) 


the  limit  can  now  be  written 


1 


lim  „ .. 

1V-.00  2 N  -i 


|  £  e-i2»r/iT0 
'fc=— // 


(3.59) 


However,  the  expression 


sin  (2A^  +  l)ff/r0 
sin  7t/T0 

is  periodic,  with  period  T0.  So,  the  limit  of  the  above  summation  becomes 


(3. GO) 


1 


lim  . 

N-+co  2N  + 1 


JV 

E* 

k=—N 


—j2xfkT0 


'mr£'('-r)  (3M) 

J°  n=0  V 


Substituting  the  above  expression  into  Eq  (  3.45)  the  average  power  spectral  density  becomes 


Sy(f)  =  yo  1^(«2)  -  £2(«)  I'M/)!2  +  Y0E^a) l$£(/)|2 £ 6  {f  ~  i)  }  (3,62) 

This  expression  for  the  power  spectral  density  is  similar  to  those  derived  by  Lawson  and 
Uhlenbeck  (11:43-44).  However,  this  expression  takes  into  account  both  random  amplitudes  and 
random  intervals,  where  Lawson  and  Uhlenbeck  have  separate  equations  for  each  case.  Also,  the 


3-16 


magnitude  of  the  expression  derived  here  is  one  half  the  magnitude  of  the  expressions  derived  by 
Lawson  and  Uhlenbeck.  The  difference  is  because  Lawson  and  Uhlenbeck  use  a  one-sided  spectrum, 
whereas  the  spectrum  derived  here  is  two-sided. 

3.2.2  Special  Cases  In  this  section,  two  special  casesof  PSD  of  a  series  of  pulses  with  random 
amplitudes  and  spacing.  First,  the  power  spectral  density  of  the  random  binary  transmission  is 
derived.  This  derivation  is  followed  by  the  development  of  the  boxcar  spectrum. 

3.2.2. 1  Random  Binary  Transmission  One  special  case  of  the  expression  for  the  power 
spectral  density  of  pulses  with  random  amplitudes  and  spacing  is  the  random  binary  transmission 
as  defined  by  Papoulis  (14:294,  341).  The  amplitudes  are  equally  likely,  and  distributed  according 
to 


P(a  =  l)  =  P(a  =  -l)=i 
Which  leads  to  the  following  expected  values: 


(3.63) 


E(a)  =  0 
E(a2)  =  1 

The  Fourier  transform  of  a  single  square  pulse  is  given  by 


(3.64) 


G(/)  =  T 


sin  (7 r/T) 

Tff 


(3.65) 


where  T  is  the  pulse  duration.  In  this  case  T  is  a  fixed  value.  Therefore,  the  deviation,  e  is  zero. 
This  implies  that  p(e)  -  6(e),  and  <!><(/)  is  unity.  From  Eq  (  3.62)  this  leads  to 
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(3.66) 


S(f)  =  £  |G(/)|2 


1  00 

^(a2)-^2(a)  +  -£;2(a)^^(/-^) 

n=0 


Substituting  Eqs  (  3.64)  and  (  3.65)  into  the  above  equation  yields 


sin  Qr/r)l2 
*7  . 

sin  2(tt/T) 
tt2/2T 


(3.67) 


This  equation  for  the  power  spectral  density  is  illustrated  in  Figure  3.3.  This  expression  agrees 
with  the  power  spectral  density  for  the  random  binary  transmission  presented  by  Papoulis  (14:341). 


Figure  3.3.  Random  Binary  Transmission  PSD 


3. 2.2. 2  Boxcar  Spectrum  The  special  case  of  the  boxcar  spectrum  is  studied  here.  The 
spectrum  considered  here  is  the  spectrum  of  the  output  of  a  boxcar  generator  (more  commonly 
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referred  to  as  a  sample-and-hold  circuit).  The  output  of  the  boxcar  generator  is  a  series  of  square 
pulses  with  random  amplitudes,  distributed  according  to  p(a),  and  with  constant  width  T .  The 
Fourier  transform  of  a  single,  unit  amplitude  square  pulse  of  width  T  is 


As  was  the  case  for  the  random  binary  transmission,  the  pulse  width  is  a  fixed  value,  and  the 
deviation,  e,  is  distributed  according  to 


p(c)  =  S(e)  (3.69) 

which  implies  that  $£(/)  =  1.  Substituting  this  value  for  $£(/)  and  Eq  (  3.68)  into  Eq  (  3.62)  the 
power  spectral  density  becomes 


S(f)  =  T 


The  following  relationship  exists  for  various  values  of 


(3.70) 


sin  7T  JT 
7T  fT 


< 


1 

0 


n  =  0 

n  =  1,2,3,  ••• 


(3.71) 


The  above  relationship  implies  that  the  delta  function  in  the  equation  for  the  PSD  only  applies 
when  /  =  0.  Therefore  the  PSD  is  written 


S(f)  =  { |£(«2)  -  £!(0)J  +  £2(<.)i«(/)| 


(3.72) 
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where  E(a2)  -  E2(a)  is  recognized  as  the  variance  of  the  amplitude  a.  This  spectrum  agrees  with 
the  power  spectrum  for  this  case  as  presented  in  Lawson  and  Uhlenbeck  (11:274). 

S.S  Pulses  with  Random  Amplitudes ,  Durations,  and  Phases 

In  the  previous  section,  the  power  spectral  density  of  a  series  of  pulses  with  random  amplitudes 
and  durations  was  developed.  In  this  section,  the  power  spectral  density  for  the  more  general  case 
of  a  series  of  pulses  with  random  amplitudes,  durations  and  phases  is  developed.  The  development 
here  follows  closely  that  of  Lawson  and  Uhlenbeck  (11:44-46).  A  single  pulse  is  defined  by  the 
expression  ane*2*J0i.  Further,  the  phase  will  change  for  each  pulse,  with  p(a)  the  probability 
density  function  of  the  phase,  a*.  Therefore,  at  the  following  times  the  pulse  y(t)  is  defined  by 


ty  :  y(t)  =  a1ej27tV°Hai') 
t2  :  y(t)  =  a2ej2*V*t+Q'+c*>') 


tn  :  7/(t)  =  aneJ’2T^+ai+a2+-+a") 


(3.73) 


where  the  au  s  are  real,  statistically  independent  random  variables  distributed  with  PDF  p(a).  The 
amplitudes  have  the  following  joint  expectations: 


E(akai)  =  < 


E(a 2) 
E\a) 


k  =  l 

k^l 


(3.74) 


The  length  of  each  pulse  is  distributed  according  to  the  PDF  p(l)  and  is  defined  by 


Ik  =  the  spacing  between  times  tk  and  4+i 


(3.75) 
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so  that  the  time  of  occurrence  of  the  kth  pulse  is 


h  =  +  h  +  h  4 - b  h-x  (3.76) 


An  example  of  a  series  of  pulses  with  random  amplitudes,  phases  and  spacing  is  illustrated 
in  Figure  3.4. 


Figure  3.4.  Series  of  Pulses  with  Random  Amplitudes,  Phases  and  Spacing 

Some  useful  relationships  presented  by  Lawson  and  Uhlenbeck  are  expanded  here  (11:45): 

©  =  ['  p(a)e>2*ada  =  A  4-  jB  =  E{e^)  (3.77) 

J~  3 
TOO 

O  =  /  dl  r-  0(/)  +  f)  =  E(e-ju'1)  (3.78) 

Jo 

1  =  r  p(l)ldl  =  E(l)  (3.79) 

Jo 
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where 


f3n  —  2iro(n 

0JC  =  2tt(/  -  /„) 


Some  relationships  involving  ft  are  presented  here.  Using  Eq  (  3.76) 


E (e>“'tk) 


eJWc<i  jjj  |eJ‘*’c0i+la4 — 


ei"ctlQk-l 


(3.80) 


since  the  /„ ’s  are  statistically  independent.  Further,  the  expectation  of  the  conjugate  of  the  previous 
equation  is 


E(e~jUctk)  =  E[(eiUctk)*] 

=  e*"e,1(fi*)*_1 


(3.81) 


Finally,  combining  Eqs  (  3.80)  and  (  3.81)  the  following  relationship  results: 


E  _ 


1  ,  k  =  l 

(ft*)'"*  ,  k  <  l 
ft*-'  ,  k>i 


(3.82) 


Likewise,  some  relationships  involving  ©  are  presented  here: 
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E  Jej(/?i+/3j+-+/?*)J  _  0k  (3  83) 

since  the  c*„’s  and  therefore  the  finys  are  statistically  independent.  Now  combining  Eq  (  3.83)  with 
it’s  conjugate,  the  following  relationship  is  arrived  at: 


1 

,  k  =  l 

E  =  <  (0*),-i 

,  k  <  l 

(3.84) 

ek-> 

< 

,  k  >  l 

A  single  pulse  from  time  tk  to  time  4+ 1  is  expressed  as 


y(t)  =  ajtei2)r(/o<+tti+«3+-+aic)  (3.85) 


The  Fourier  transform  of  this  single  pulse  from  tk  to  tfc+i  is 


Y(f]tk;tk+ i)  =  akej^l+Pl+"'+M  f  tr&'V-Wdt 

Jtk 

-  — :L.akej(Pi+Pi+-+Pi<)(e-ivctk+i  _  e-jwcU)  (3.86) 

iwc 


A  finite  series  of  N  pulses  from  time  i\  to  is 


N 

yN(t)  =  ^  ake}2*(fot+ai+a,+-+c'k)  (3.87) 

k=l 


and  the  Fourier  transform  of  this  series  of  pulses  is 
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(3.88) 


N 


YN(f)  =  - j-  ajbCi^1+^+,“+^)(e"^c<,t+1  -  e“iWctfc) 


N 


=  -rE* 

J  °  Jt  =  l 


where 

=  ajteJ(^i+P2+*--+/?fc)(e“iwcU+i  _  e"^Wc<fc)  (3.89) 

5.5.i  Power  Spectral  Density  The  average  power  spectral  density  of  the  series  is  developed 
here.  The  average  PSD  is  defined  by 


Sy(f)  =  \m^X\YN(f)\2  (3.90) 

where  Nl  is  the  average  length  of  N  pulses. 

The  square  of  the  magnitude  of  the  Fourier  transform  is 


IW)I2 


N  N 


^EE*'4? 


c  k=l  1= 1 


1  N  N-lN-k 

~2  Y2  l^l2  +  Y2,  Y2  {AWAk  +  Ak+lAk) 

c  Ln=l  k= 1  /=1 

1  /  N  N-lN-k  \ 

c  \n= 1  k-1  1=1  / 


(3.91) 


since  Z  +  Z*  =  2 Re(Z).  Using  Eq  (  3.89)  the  first  summation  in  the  above  equation  is 
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n  =  1 


^  |/ln|2  =  |  (e~-,a’c<n+l  -  e"','a’c<n)] 

n=l 

X  ^a„e_;^1+^3+”'+^")  (e;Wctn+1  -e;a’c‘n)J| 

N 

=  El°"  (e-;wc<n+i  _  e-iw«‘»)  an  (eJ‘u'c‘"+1  -  ei"‘‘»)] 
n=l 

N  r  i 

=  a\  2  —  eJ’Wc(*B+l“*n)  -  £-*"«(*«+!-*») 


(3.92) 


n= 1 


The  expected  value  of  this  summation  is 


■(!H 


^  #  ja2  ^2  -  -  e-M=(<>.+i-<>>)j  j.  (3  93) 

n=l 
N 

a2)  ^2-  E 


pwc((n+i-i„)'  _  E  Je-ia>c(«n+i-<„)|  |  (3.94) 


Using  Eq  (  3.82)  this  expectation  becomes 


N 


N 


E  £ki2  E^(a2)(2-^-fi*) 


\n  =  l 


n=l 


(3.95) 


Since  Q,  4-  ft*  =  2 i?e(ft) 


*(JH 


NE(a2)  [2  -  2iZe(fi)] 
2N E(a2)  [1  —  <j>(f)] 


(3.96) 


Using  Eq  (  3.89)  the  terms  inside  the  second  summation  are 
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Ak+iA*k  =  [aHi^^,+^+"'+^+,)(c"jWeU+,+1 -e_iW4U+')] 

x  ■+0fc)(giwcU+i  _  cjwctfc)j 

=  ak+iakejWk+l+Pk+>+'''+Pk+l) 

x  ^e;w,^k"“^+I)  —  eiwc(^“^+*+0  —  c/w«(<*+i -<*+*)  ejwc(4+i-U4i*n)J  (3.97) 
The  expected  value  of  this  term  is 


J?(4k+Ml)  =  £(a*+,a*)£  [e#tf»+i-Ww+-+A+i)]  {j?  j^c(<fc-u+,)j 

_£  |eJw«(,k-tk+l+»)j  _£  |g»«c(<«f+l-U+l)j  +  £  |fiJ«e(tl(+»-<*+l*l)j  J.  (3.98) 


By  using  Eqs  (  3.74),  (  3.82),  and  (  3.84)  this  expectation  can  be  expressed 


E(Ak+iAi)  =  jEB(a)e,(n,-n,+1-n,”1+n') 

E{Ak+iAl)  =  E2{a)  2 (0©)' -  0(fi©)' -  i(fi©)' 


(3.99) 

(3.100) 


The  summation  of  this  expectation  is 


(iV-liV-l;  \  iV—  1  N—k  r  - 

2  E  4t+«4E  =  £2(«)  E  E  2(Q0)'  -  Q(Q0)'  -  o(fi0)' 

k= 1  i=l  /  *=1  1=1  *■ 


(3.101) 


Let  y  =  (f2G).  Now,  the  summation  can  be  written 
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N-k 

N-k-l 

EV 

=  E 

j=i 

m=0 

N-k-l 

II 

M 

5 

m= 0 

•a 

1 

* 

<*- 

1 

i— < 

II 

1-7 


(3.102) 


The  expectation  of  this  summation  now  becomes 


*(££***) = 


=  E*(a) 


7  / 

2- Q- 

h)] 

11  —  7  V 

7  / 

2  — a- 

-)} 

L  i — 7  \ 

n)\ 

N-l 


=  i?2(a)r^©("fi2+2Q'1) 

=  -^2(«)r=^©  (!  -  f2)2  -  r^-  (!  -  T^"1)]  (3.103) 


Substituting  Eqs  (  3.96)  and  (  3.103)  into  Eq  (  3.91)  the  square  of  the  magnitude  of  the 
Fourier  transform  becomes 


r 

IW)|2  =  ^|{2  NE(a2)[l-m) 

-2£S<0)Ke{rls(1-n)! 


N  - 1 


1-7 


(l-T*”1) 


(3.104) 
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Using  Eq  (  3.90)  the  power  spectral  density  is  obtained: 


Sy(f)  =  2h2{f  _—oji  j^0’)  t1 "  WM  "  ^  ^Re  [l^Q  tt  -  fi)2]  }  (3.105) 


The  above  equation  for  the  power  spectral  density  can  be  quite  useful  in  the  given  form, 
especially  when  the  terms  ©  and  ft  are  real  quantities.  However,  this  equation  can  be  expanded  as 
a  function  of  the  component  parts  of  ©  and  ft.  The  coefficient  of  the  E2(a)  term  can  be  written 


Re 


© 

i-fte 


(1-fi)2 


Re  [©  (1  -  ft)2  (1  -  ft©)’ 

~  |l-ft©|2 

N_ 

D 


(3.106) 


where 

N  =  the  numerator 
D  =  the  denominator 

Using  Eqs  (  3.77)  and  (  3.78),  the  denominator  in  the  above  equation,  D)  can  be  written 


D  =  |l-ft©|2 

=  (l-ft©)(l-ft©)* 

=  {[1  -  M{f)  +  Di>(f)\  -  j  [Aip(f)  +  B<t>{!))} 
x  {[1  -  AW)  +  +  j  [AW)  +  **(/)]} 

=  [1  -A<Kf)  +  rhp(f))2  +  [Mf)  +  B<Kf)}2  (3.107) 
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The  components  of  the  numerator  are  expanded  here: 


(i-ft)2  =  [i 

=  [1  -  2 4(f) + fa)  -  -  mi)  [i  -  m 

0(1  -  ft)2  =  (A  +  jB )  { [1  -  2 4>(f)  +  <j>2(f)  -  i>2(f))  -  j2xp(f)  [1  -  «>(/)]} 

=  {A  [1  -  2<p(f)  +  <&2(/)  -  V-2(/)]  +  2 BiP(f)  [1  -  } 

+j  {B  [1  ~-2m  +  <f,2(f)  -  tfifj]  -  2  AW)  [1  -  <f>(f))}  (3.108) 

(1-ft©)*  =  {l +  jrp(f))(A  +  jB)y 

=  [1  -  A<f>(f)  +  BxP(f))  +  j  [Ai>(f)  +  B<f>(f))  (3.109) 


Substituting  Eqs  (  3.108)  and  (  3.109),  the  expression  for  the  numerator,  N ,  can  be  written 


N  =  Re  (1  —  ft)2  (1  —  ft©)*] 

=  [1  -  A<j>(f)  +  B'ip(f))  {A  [1  -  2(j>{f)  +  <f>2(f)  -  i>2{f))  +  2 BiP(f)  [1  -  <f>(f))} 

- [Ai>(f)  +  Bt(f)}  {B  [1  -  2 +  <t>V)  ~  V>2(/)]  -  W{f)  [1  -  <£(/)]} 

=  [1  -  2 0(/)  +  t2(f )  -  i>2(f)}  [A  -  A2<f>(f)  -  B2<f>(f)] 

+  [1  -  *(/)]  [mU)  +  2 B2i>2(f)  +  2A2tf(f)} 

=  [1  -  2 +  4>\f)  -  {>1  -  <!>(/)  [A2  +  B2} } 

+2 V-(/)  [1  -  m\  {B  +  m  [A2  +  B 2] }  (3.110) 


Substituting  Eqs  (  3.107)  and  (  3.110)  into  Eq  (  3.106),  and  after  some  simplification,  the 
power  spectral  density  becomes 
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Sy«)=2 

(1-2^+^-^) [yi-^(>i2+^2)]+2^(i-^)[B+^(/i2+^2)n  rt1in 

(1  -A<f>  +  B^)2  +  (A^  +  B<f>)2  )  {' 

where  the  <£  and  rp  terms  are  in  general  functions  of  frequency,  /.  The  above  equation  does  not 
reduce  as  elegantly  as  the  equation  presented  in  Lawson  and  Uhlenbeck  (11:45).  This  is  due  to 
the  introduction  here  of  random  amplitudes,  while  Lawson  and  Uhlenbeck  only  consider  random 
phases  and  pulse  spacing. 

3.8.2  Special  Cases  The  power  spectral  density  presented  can  be  a  very  useful  tool  to  derive 
the  power  spectra  of  various  waveforms.  Two  special  cases  of  random  waveforms  are  considered 
here.  First,  the  PSD  for  the  random  binary  transmission  is  derived.  Finally,  the  power  spectra  for 
the  random  telegraph  signal  is  developed. 

3.3.2. 1  Random  Binary  Transmission  The  special  case  of  power  spectral  density  of 
the  random  binary  transmission  is  derived  here  from  the  general  case  of  the  PSD  of  a  series  of 
pulses  with  random  amplitudes,  spacing,  and  phases.  Since  the  random  binary  transmission  is  a 
baseband  process,  f0  is  set  to  zero.  Two  methods  to  achieve  this  PSD  are  illustrated  here. 

First,  the  amplitudes  are  assumed  to  have  equally  likely  distributions  as  follows: 

P(a=l)  =  P(a  =  -l)=i  (3.112) 

which  is  the  same  as 


?(«)  =  \  P(«  + 1)  +  <5(a  -  1)] 


(3.113) 
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so  that  the  expected  values  are  given  by 


E(a)  =  0  (3.114) 

E{a2)  =  1  (3.115) 


The  phase  will  be  assumed  to  be  a  constant  value  of  zero.  This  is  because  the  random  amplitudes 
account  for  a  phase  change  of  either  0  or  n  radians  (with  an  equally  likely  distribution).  Substituting 
the  above  expectations  int  Eq  (  3.105) 


S(f)  = 


[W(/)) 

2br2/2 


(3.116) 


For  the  random  binary  transmission,  the  pulse  spacing  is  a  constant  value,  T,  so  that 


p(l)  =  S(l-T) 


(3.117) 


and  7  =  T.  Next  <f>(f)  is  derived  from  p(l)  by  use  of  Eq  (  3.78)  as  follows: 


Q 


/‘OO 

/  p{l)e-j2n,f  dl 

Jo 

rOO 

/  p(l)  =  8(l-T)e~j2*,fdl 

Jo 

e-j2*JT 


=  cos  2tt  fT  +  j  sin  2tt  fT 


(3.118) 
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Now,  =  cos  2 tt/T,  and  the  PSD  is  now 


S(f) 


(1  —  cos  2tt  JT) 

m&p 

2  sin2  vfT 

22V/2 

/sin7r/r\2 

V  *fT  ) 


(3.119) 


which  is  of  the  same  form  as  the  PSD  of  the  binary  random  transmission  presented  by  Papoulis 
(14:341).  This  power  spectral  density  is  illustrated  in  Figure  3.3. 

A  slightly  different  method  will  be  demonstrated,  which  yields  the  same  result.  In  this  case, 
take  the  amplitudes  to  be  a  constant  value,  so  that 


E(a )  ss  E(a2)  =  1 


(3.120) 


and  let  the  phase  changes  be  either  0  or  it  radians,  with  equal  likelihood,  so  that 


<5(a)  +  6 


From  Eq  (  3.77)  the  value  of  ©  is  derived: 


(3.121) 


0 


p(cv)e;27rorda' 


/_■  i  [«<.)+ 


0 


(3.122) 
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and  A  =  B  =  0,  so  both  Eq  (  3.105)  and  Eq  (  3.111)  reduce  to 


*(/)  = 


[i-m 

2/tt2/2 


(3.123) 


Since  the  probability  density  of  the  pulse  spacing,  p(/),  has  not  changed,  the  PSD  reduces  to  the 
form  of  Eq  (  3.119). 


3.S.2.2  Random  Telegraph  Signal  In  this  section  the  PSD  of  the  random  telegraph 
signal  is  derived.  This  waveform  is  defined  by  Papoulis  as  a  point  process  with  an  underlying 
Poisson  distribution  (14:288-290).  This  process  is  a  baseband  process,  so  the  frequency  f0  will  be 
set  to  zero.  The  interarrival  times  of  this  point  process  will  determine  the  pulse  length,  /,  for  the 
case  considered  here.  The  interarrival  times  of  a  Poisson  process  are  shown  by  Davenport  to  have 
an  exponential  density  (3:471).  Therefore,  in  this  case  the  probability  density,  p(/),  is  given  by 


p(/)  =  Ae“A?,/>0  (3.124) 

For  the  random  telegraph  signal,  the  pulse  amplitudes  transition  from  -1-1  to  -1,  and  vice 
versa,  at  each  interarrival  time.  Therefore,  one  way  to  view  the  distribution  of  the  amplitudes  and 
phases  is  to  set 


E{a)  =  E(a2)  =:  1 


(3.125) 


and  the  density  of  a  to 


p(cv)  -  5 


From  Eq  (  3.77)  this  leads  to 


(3.126) 
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©  =  A  +  jB 


(3.127) 


so  that  1  =  —1  and  B  =  0.  Substitution  of  the  above  expressions  into  Eq  (  3.111)  yields  the  power 
spectral  density  of  the  signal: 


S(f)  = 


1 


2/7r2/2 

1 

h2P 


(W)- 
l-p-p 


(1-2 ?S2  -  ^2)  (-1  -  <ft)  +  2^2  (1  -  <f>) 

(l  +  tf  +  P 


(l  +  <f>)  +  p 


(3.128) 


From  Eq  (  3.124),  the  following  expectation  results: 


r co 

7  =  /  l\e~xldl 

Jo 

1 

~  \ 


(3.129) 


And,  from  Eq  (  3.78) 


ft  =  + 

=  /%( l)e~i2*J,dl 

Jo 
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/»00 

ft  =  /  \e~lV2*J+xdl 

Jo 

A 

A  +  jw 


where  «  =  2~ /.  ft  can  be  broken  onto  its  component  parts  as  follows: 


ft  = 


A2  .  Aw 
A2  +  w2  JA2  +  w2 


Substituting  the  above  into  Eq  (  3.128),  the  PSD  is  obtained: 


(3.130) 


(3.131) 


S(w) 


S(f) 


4A 

1  X4  AV 

A  (A3+w3)3  (A3+w3)4 

u>2 

Ui.  \2 

4A 

L  V"  ‘  A3+w3  J  '  (A3+w3)' 
(A2  +w2)2  -  A4  -  A2w2 

w2 

(2A2  +  u>2)2  -f  A2w2 

4A 

■  w4  +  A2w2  • 

w2 

4A4  4*5A2o;2  + 

4A 


4A2  +  w2 
4A 

4A2  +  (2tt/)2 


(3.132) 


This  equation  for  the  power  spectral  density  is  illustrated  in  Figure  3.5. 

This  expression  is  the  same  as  that  given  by  Papoulis  for  the  power  spectral  density  of  the 
random  telegraph  signal  (14:341). 
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Figure  3.5.  Random  Telegraph  Signal  PSD 
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IV.  Probability  Densities  and  Expectations 


The  probability  density  functions  of  the  noise  waveform  and  the  signal  plus  noise  waveform, 
along  with  the  power  spectral  density,  are  one  of  the  most  important  properties  required  in  the 
analysis  of  the  effects  of  noise  jamming  on  radar  receivers.  The  PDF  describes  the  distribution  of 
the  amplitude  of  the  noise  jamming  waveform.  In  this  chapter  the  probability  density  functions  of 
the  amplitude  of  the  output  of  both  a  linear  and  quadratic  detector  are  derived.  These  PDF’s  are 
derived  for  the  case  of  a  noise  only  input  as  well  as  the  case  of  the  input  of  signal  plus  noise.  Lastly, 
several  expected  values  are  derived.  The  mean  and  mean  square  values  are  determined  for  the  two 
cases:  noise  only  and  signal  plus  noise.  Also,  the  autocorrelation  function  for  the  signal  plus  noise 
case  is  derived  for  both  detector  types. 

f.l  Probability  Densities  of  the  Output  of  a  Linear  Detector 

In  this  section,  the  probability  density  functions  for  the  output  of  a  linear  (envelope)  detector 
are  developed.  The  PDF’s  will  be  derived  for  the  case  of  the  input  of  noise  only  as  well  as  the  case 
of  a  signal  plus  noise. 

Probability  Density  for  Noise  Only  The  noise  analyzed  here  is  assumed  to  be  white 
and  have  a  Gaussian  distribution  with  mean  value  of  zero  and  a  variance  of  a2.  This  noise,  n,(t), 
is  input  to  the  IF  filter,  whose  transfer  function  is  H(f).  The  noise  output  from  the  IF  filter,  n(2), 
can  be  portrayed  in  quadrature  form  as 


n(t)  =  nx(t)  cos27rfct  +  nfJ(t) sin2nfct 


(4.1) 


where 


fc  =  the  center  frequency  of  the  IF  filter  passband 
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The  noise  components  nx(t )  and  ny  (t)  are  independent,  zero  mean  Gaussian  random  variables  with 
PDF’s 


P(nx) 

pK) 


i 


V2tcW 

1 

y/2 tiW 


exp  - 


exp  - 


2  W 

2  W 


(4.2) 

(4.3) 


where 


/OO 

\H(f)\ 2  df 

•CO 


(4.4) 


Because  these  two  noise  components  are  statistically  independent,  their  joint  PDF,  p(nx,ny),  is 


,  X  1  f  *1+"*' 


(4.5) 


The  output  of  the  linear  detector,  L(t ),  is  the  envelope  of  the  input  signal,  defined  by  the 
expression 


Ln(t)  =  y/nl(t)  +  nl(t)  (4.6) 

where  the  subscript  n  is  to  denote  the  noise  only  condition.  The  terms  nx  and  nv  can  be  transformed 
to  polar  form  in  terms  of  Ln  and  <j)n  where 


nx  —  Ln  cos  (f>n  — 

(4.7) 

fly  =  Lfx  sin  <f)n  =  1^2  $n) 

(4.8) 
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The  Jacobian  of  this  transformation  is 


dh  i 

Qhy 

OLn 

d<f>n 

dhx 

Qhn 

dLn 

d$n 

which  reduces  to 


(4.9) 


J  =  Ln  cos2  <j>n  +  Ln  sin2  <f>n 
=  Ln 


(4.10) 


The  joint  density  of  Ln  and  <j>n  can  be  derived  from  the  joint  density  of  nx  and  ny  by  the  equation 


P(Ln,<f>n) 


Pn„n v  [hi  (L  n  >  $n)  >  ^2  (I'm  $n)]  J 

1  f  Ll  cos2  <f>n  +  Ll  sin2  <j>n  \ 
2 nW  GXP  V  2 W  J 


Ln 


(4.11) 


From  the  joint  PDF,  the  marginal  PDF  of  the  envelope,  p(Ln) ,  can  be  obtained  by  integrating  over 
all  values  of  <j>n  so  that 


ff 

0 


(4.12) 
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for  all  values  of  Ln  >  0.  This  probability  density  is  known  as  the  Rayleigh  PDF  and  is  illustrated 
in  Figure  4.1. 


Figure  4.1.  PDF  of  Noise  Only  for  a  Linear  Detector 


It  can  be  shown  in  a  similar  manner  that  the  phase  term,  <f>U)  is  uniformly  distributed  over 
the  interval  (0,27t),  and  that  Ln  and  <j>n  are  statistically  independent  random  variables. 

Jf.1.2  Probability  Density  for  Signal  Plus  Noise  The  PDF  for  a  signal  plus  Gaussian  noise 
will  now  be  derived  for  the  case  of  a  linear  detector.  It  is  assumed  that  the  signal  and  noise  are 
added  together  prior  to  the  detector,  and  that  the  signal  is 


s(t)  =  a(t)  cos  2tt/c*  +  f3(t)  sin  2t: fct  (4.13) 

so  that  from  Eqs  (  4.1)  and  (  4.13)  the  input  to  the  detector,  e(t),  is 
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e(t)  =  s(t)  +  n(t) 

=  [a(t)  +  nx(t)}  cos  2irfct  +  \fl(t)  +  nv(t)]  sin  2i xfct 

=  x(t)  cos2tt  fct  +  y(t)  s\n27r  fct  (4.14) 

where  a:(t)  and  y(t)  are  independent  Gaussian  random  variables.  Since  nx(t)  and  %(<)  are  zero 
mean  random  variables,  it  follows  that  the  means  of  x(t)  and  y(t)  are  now 


E[x(t)] 

=  «(*) 

(4.15) 

E[y(t)) 

=  m 

(4.16) 

and  the  variance  of  each  is  W ,  as  was  the  case  for  noise  only.  Now,  the  joint  PDF  of  x  and  y  can 
be  written: 


(s  -  Of)2  +  (y  -  /?)2 


2  W 


(4.17) 


This  density  can  be  transformed  into  polar  form  in  terms  of  Ls+n  and  <^+n,  as  was  done  for  the 
case  of  noise  only,  using  the  equation 


PiLs+n  j  fis+n)  —  Px,y  [^1  ^2 >  ^i+n)]  J  (4. IS) 
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where 

^l(^i+n )  =  -kj+n  COS^+n 

>  ^5+n)  =  ^s+nsi^^3+n 

J  =  Lj+n 

This  leads  to  the  joint  density 


p(Ls+n,<f>a+n)  -  |j|£exp 


T'a+n 

27TPF 


exp 


(L,+n  COS  <ft,+n  ~  «)~  +  (7;,+n  sin  <ft,+n  -  /?)2 
217 

J-%±n  ~  2L±+n  (tv  COS <f>3+n  +  /?sin  ^.fn)  +  a2  +  /?2 


2W 


(4.19) 


To  arrive  at  the  PDF  of  the  amplitude  of  the  detector  output,  p(L3+„),  the  above  joint  PDF  is 
integrated  relative  to  <f>3+n  over  the  interval  (0, 2sr): 


Pi^s+n)  =  /  p(-^3+«  j  ^4+»)^i+n 

JO 


/*  2x 

Jo 


^s+n 


exp 


7^+n  -  2L,+n  (a  cos  </>,+„  +  /?  sin  0J+n)  +  a2  +  (32 1 


2ttV7 

7-i+n  /  +  <*2  +  /?2 


2VF 


•fn 


»°w(-Z*±{£±£.} 


2irW 
/» 

x  /  exp 
Jo 


■^5+n 


2V7 

(a  cos  ^+„  +/?  sin  <£,+„) 


d<f>i+n 


(4.20) 


The  integral  in  the  above  equation  is  of  the  form  of  Eq  (3.937,2)  found  in  Gradshteyn  and  Ryzhik, 
and  reduces  to  (8:488) 


Ls+n 

W 


(acos<j>3+n  +/?  sin  <£,+„) 


d<f>s+n  =  2^7,, 


(4.21) 
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and  upon  substitution  of  the  above  equation  into  Eq  (  4.20),  the  PDF  of  the  envelope  now  becomes 


K^+n) 


2nW  6XP  i" 

Lj+n  ( 
exp  \ 


+  CV2  +/?2\ 

2iy  y  0 


(4.22) 


for  >  0.  The  term  I0  in  the  above  equation  is  the  modified  Bessel  function  of  order  zero.  This 
expression  for  the  probability  density  is  known  as  the  Rician  PDF.  The  Rician  PDF  is  illustrated 
in  Figure  4.2,  similar  to  the  form  of  a  related  plot  found  in  Lawson  and  Uhlenbeck  (11:154).  Two 
values  of  z  =  (a2  +  /32)/2 W  are  plotted  in  Figure  4.2. 


Figure  4.2.  PDF  of  Signal  Plus  Noise  for  a  Linear  Detector 


DiFranco  and  Rubin  provide  useful  approximations  for  two  cases  of  the  modified  Bessel  func 


tion.  For  the  small  argument  case  (5:344) 
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(4.23) 


.  ,  .  .  x2  x* 

I0{x)  =  l  +  — +  —  +  ■“, X<1 

and  for  the  large  :  -pument  case  (5:344) 

iW“^(1  +  C-)-*>1  (1'M) 

JJ.3  Expected  Values  The  expected  values  of  L  and  L 2  for  the  noise  only  and  signal  plus 
noise  cases  are  derived  here*  Also,  the  autocorrelation  function,  L1L2  for  the  signal  plus  noise  case 
is  approximated.  First,  the  case  of  the  average  value  of  the  envelope  of  the  noise  only  output  is 
considered.  With  the  PDF  of  Ln  given  by  Eq  (  4.12),  the  average  value  of  the  noise  amplitude  is 
given  by 


/00 

LnP{Ln)dLn 

■co 

=  2 a  [  Lle~aL*dln  (4.25) 

Jo 


where 

a“  2*7 

The  above  integral  is  of  the  form  of  Eq  (666)  in  the  CRC  Standard  Mathematical  Tables  (16:465). 
Using  this  relationship,  the  expectation  becomes 

E(Lny=^Y  (4-26) 

Next,  the  expected  value  of  L„  is  calculated: 


4-8 


(4.27) 


/CO 

LlV{Ln)dLn 

=  l  w^{-m)dL’' 

~  2W  J0  ”eXP(~2 w)da 

where 

a  =  L2n 
da  =  2  LndLn 


This  integral  is  of  the  form  of  Eq  (661)  in  the  CRC  Tables  (16:465).  After  substitution,  the 
expectation  becomes 


=  2  W  (4.28) 

Now,  the  expectations  for  the  signal  plus  noise  case  are  derived.  With  the  PDF  of  the  signal 
plus  noise  output  given  by  Eq  (  4.22),  the  expectation  of  is 


Ls+nP(Ls+n)dL$+n 


E(Ls+n)  =  J° 

■  r%-( 


tj+n  +  Or2  +  /?2 
2  W 


ell 


5+n 


(4.29) 


The  integral  in  the  above  equation  is  defined  in  Lawson  and  Uhlenbeck  by  the  following  expression 
(11:174): 
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(4.30) 


r(f)  ($)‘ 

2pi‘T(i'+  1) 


exp 


«■"(=? 


—  /i  ,  „  „  a 

+  1>1/+1;"v 


For  the  case  of  u  =  0,  the  above  equation  becomes 


(4.31) 


where  iFt(a,6;«)  is  the  confluent  hypergeometric  function,  defined  by  (11:174) 


\Fi(a,b;z) 


a  a(a  +1)  a(a  +  l)(a  +  2)  z3 

1+b  +  6(6+1)  2!  +  b(b  +  l)(b  +  2)  3!  +"' 


(4.32) 


Using  the  integral  relationship  of  Eq  (  4.31),  the  expectation,  E(Ls+n),  becomes 


E{L3+n) 


V  2  W  ) 


£  (I)  (2W)§ 

2 


exp 


(<S±£_\ 

V  2W  ) 


iFi 


(  1  1  <*2  +  /?2\ 
v  2’  ’  2W  ) 


(4.33) 


From  the  CRC  Tables  Eq  (605),  T  (|)  =  ^  (16:460).  With  this  relationship,  the  expectation  can 
be  reduced  to  the  following  expression: 


F(Ls+n)  — 


Q2  +  /?2  \ 

2  W  ) 


(4.34) 


Lawson  and  Uhlenbeck  provide  an  asymptotic  approximation  for  \  F\  (a,  b;  z)  for  large  negative 
values  of  z  (11:174).  This  approximation  is 
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iFi(a,b-,z)zi 


m 

T(b  -  a) 


(-*)-“ 


a(a-b  +  1)  (  a(a  +  l)(a-6  +  l)(a-6  +  2)  , 

1 - - -  +  - "o - + 


(4.35) 


In  the  above  case  for  E(Ls+n),  where  z  = 


the  confluent  hypergeometric  function  becomes 


iFi 


\  „  r(i) 
!~  r(|) 


M 


1 


3 


1 


(4.36) 


By  keeping  only  the  first  term  of  the  expansion,  and  substituting  into  Eq  (  4.34)  the  average  value 


of  the  detector  output  is  approximated  by 


E(Ls+n)  « 


pxWT(l)  far 

V  2  r(|)  V 


+P2 


T  (|)  V  2 W 
nW  2  /a2+/?2 

2  v^V  2.1V 


«  \A2  +  /32 


(4.37) 


Now,  the  expected  value  of  L2+„  is  derived.  Using  Eq  (  4.22)  the  expectation  is  written 


/oo 

■oo 


-oo 

„  r  -r  Z'  ^+n+«2  +  /?2\  T  /%+„  r^rrA  „  „  _ 

-  J0  —  exP^  2W  ) Io  \~w~Va  +0J  dL‘+n  (4-38) 


From  Eq  (  4.31)  the  mean  squared  value  becomes 
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„,r2  v  _  1  ...J  a2  +  p2\\V(2)4 W2__fa2  +  p2\ 

E(L,+n)  VFeXPV  2 W  j[~  2  exp(  2W  )lFl\  1,1; 


a2  +  /?: 
2W 


2\1 


=  2WiFi  I  —1,1;  — 


a2  +  ft2 
2W 


(4.39) 


It  follows  from  Eq  (  4.32)  that  the  confluent  hypergeometric  function  in  the  above  equation  reduces 
to 


1^1  (—i,  l;  z)  =  I-* 


(4.40) 


Therefore,  the  mean  square  value  becomes 


E(LUn)  =  2w(l  +  ^P) 
=z  2  WF  +  <*2-f/?2 


(4.41) 


The  last  expectation  considered  for  the  linear  detector  output  is  the  autocorrelation  function 
of  the  amplitude  for  the  case  of  signal  plus  noise.  First,  the  joint  PDF  of  x\}  y\}  a; 2  and  1J2  is 
required.  The  subscripts  refer  to  the  times  and  *2  respectively.  This  PDF  is  a  four-variate 
Gaussian  of  the  general  form 


p(X)=  (2tt)S|V|5  exp 


--(X-mfV-l(X-m) 


(4.42) 
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where 


■ 

* 

XI 

<*i 

Vi 

;  m  = 

h 

<*2 

V2 

J2  . 

V  is  the  covariance  matrix  defined  by 


1  0  p  0 


V  =  W 


0  10  p 

p  0  10 


0  p  0  1 


(4.43) 


which  follows  from  the  fact  that  the  X{  and  y,*  terms  are  uncorrelated,  but  x\  and  X2  are  not 
uncorrelated.  Likewise  for  y\  and  ?/2*  The  term,  p,  is  the  normalized  correlation  coefficient  defined 
by 


__  cov(a;,',a;j) 
P  ~  a(xi)cx(xj) 


(4.44) 


As  Lawson  and  Uhlenbeck  point  out,  the  equation  for  the  normalized  correlation  coefficient  in  this 
case  can  be  written  in  terms  of  the  IF  filter  transfer  function,  H(f),  as  follows  (ll:155)f 


,  x  fZoUKf)?  ^os2~  frdf 

P[T)  ~  IZ 


where 


T  =  t2-tl 


(4.45) 
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The  above  leads  to  the  following  equation  for  the  joint  PDF  of  y*,  #2  and  y2‘ 


_  P2)  P2)  {(a:i  ~  Ql)2  +  (yi  ~  A)2  +  (*2  ~  az)2 


4;pW2(1 

+(j/2  -  /?2)2  -  2/>[(xi  -  ai)(a;2  -  <*2)  +  (yi  -  /?i)(?/2  -  fa)]  j 


(4.46) 


Now  the  autocorrelation  function  can  be  derived: 


L,L,  =  E  +  v'i'j  + 

=  +  »/?)  (a/^+I/2)  P{xi,yi\X2,y2)dxldyldx2dy2  (4.47) 


The  integrals  in  the  above  equation  can  not  be  solved  exactly.  However,  Lawson  and  Uhlenbeck 
provide  two  approximations  for  these  integrals:  one  for  the  small  signal-to-noise  ratio  case,  and 
one  for  the  large  signal-to-noise  ratio  case  (11:156-157).  For  the  small  signal-to-noise  ratio  case  the 
autocorrelation  is  approximately 


LlLi  *  3^7  [K  +  $  +  4W)  («j  +  Pi  +  W) 

-\-AW  p{ot\ai  +  /?i#>)  +  4W2p2]  (4.48) 


and  for  a  large  signal-to-noise  ratio  the  autocorrelation  is  approximated  by 
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L\L>2  ~  [(<*1  +/?l)  (<*2  +/??)] **  +  Wp 
,  (aia2+A/?2)2 


ai»2  +  Pi  P2 


[M  +  fl)(«I+tf)P 


[(a?  +  /??)(^+/?22)]! 


(4.49) 


Probability  Densities  of  the  Output  of  a  Quadratic  Detector 

The  PDF’s  of  the  output  of  a  quadratic  (square-law)  detector  are  developed  in  this  section. 
As  was  the  case  for  the  linear  detector,  the  PDF’s  for  both  noise  only  and  signal  plus  noise  are 
presented  here.  The  output  of  the  quadratic  detector,  Q(i),  can  be  thought  of  as  merely  the  square 
of  the  output  of  the  linear  detector.  For  this  case,  Papoulis  gives  an  expression  for  the  PDF  of  the 
output  of  a  square-law  detector  for  y  =  ax2: 


where 

px(x)  =s  the  PDF  of  the  random  variable  x 
U(y)  =  the  unit  step  function 


(4.50) 


Since  the  PDF’s  p(Ln)  and  p(L,+n)  are  non  zero  only  for  L  >  0,  the  PDF  of  the  output  of  the 
quadratic  detector,  Q(t)  =  L2(t),  can  be  written 


*q)ss23qp*W*)  (4-51) 

This  result  will  be  used  to  derive  the  PDF’s  for  the  quadratic  detector.  As  before  the  signal  out  of 
the  IF  amplifier  is  given  by  Eq  (  4.13)  and  the  noise  process  is  given  by  Eq  (  4.1). 
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For  the  case  of  noise  only,  the  PDF  for  the  linear  detector  was  given  by  Eq  (  4.12).  After 
applying  Eq  (  4.51),  the  PDF  for  the  output  of  the  quadratic  detector  for  an  input  of  Gaussian 
noise  only  is 


(4.52) 


which  is  only  non  zero  for  Qn  >  0.  Therefore,  the  output  of  the  quadratic  detector  has  an  expo¬ 
nential  probability  density.  This  PDF  is  illustrated  by  Figure  4.3. 


Figure  4.3.  PDF  of  Noise  Only  for  the  Quadratic  Detector 


A  similar  approach  is  taken  for  the  signal  plus  noise  case.  By  applying  Eq  (  4.51)  to  the  PDF 
of  the  signal  plus  noise  given  by  Eq  (  4.22)  the  PDF  of  the  signal  plus  noise  output  of  the  quadratic 


4-16 


detector  p(Q5+n)  is  obtained: 


p{Qs+n) 


Figure  4.4  illustrates  this  probability  density  function  for  two  separate  values  of  z  =  (a2+/?2)/2W/, 
similar  to  to  plot  of  the  PDF  for  the  case  of  the  linear  detector,  Figure  4.2. 


Figure  4.4.  PDF  of  Signal  Plus  Noise  for  a  Quadratic  Detector 


4.2.1  Expected  Values  The  expected  values  for  the  output  of  a  quadratic  detector  are  pre¬ 
sented  in  this  section.  The  expectations  of  Q  and  Q 2  for  the  cases  of  noise  only  and  signal  plus 
noise  are  derived.  The  autocorrelation  function  for  the  signal  plus  noise  case,  is  also  derived. 
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The  expectations  for  the  noise  only  case  are  considered  first.  With  the  probability  density  of  Qn 
given  by  Eq  (  4.52)  the  average  value  of  the  quadratic  detector  output  with  only  a  noise  input  is 


E(Qn) 


Qnp{Qn)dQn 


dQn 


(4.54) 


The  integral  in  the  above  equation  is  of  the  form  of  Eq  (661)  from  the  CRC  Tables  (16:465).  Using 
this  equation  for  the  integral,  the  expectation  becomes 


E{Qn)  =  2^(2W)2 

=  21 V 


(4.55) 


Next,  the  mean  square  value  of  Qn,  E(Q%),  is  calculated  from  the  PDF  given  by  Eq  (  4.52): 


ml) 


Qlv{Qn)dQn 

oo 


Again  using  CRC  Eq  (661),  the  mean  square  value  is  now  written 


(4.56) 


E(Ql)  =  ^(2  Wf 

=  8  W2 


(4.57) 
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Now  the  ease  for  signal  plus  noise  is  considered.  The  PDF  of  the  amplitude  of  the  output  of 
the  quadratic  detector  when  the  input  is  signal  plus  noise,  p(Qj+n)>  is  given  by  Eq  (  4.53).  This 
leads  to  the  expected  value  of  Q4+n: 


$(Qs+n) 


=  1. 


-oo 

00 


Q  s+nP{Q  $+ti)dQ  s+n 

Qa+n  f  Qs+n  +  +  P2 


r  ( 

=  L  ™rexp{- 


2W 


Io 


^/+,>  V“2+/?2  )  dQs+n 


w 


=  w^ir^w-)rt3‘xp(~^)h{^'/^TW)d‘  (“'58) 


where 

t  =  \/Q7+n" 

dt  -  vfc:dQ,+n 

From  Eq  (  4.31)  of  Section  4.1.3  the  expected  value  becomes 


^(Qs+n, 


r(2)(2  wf 

— -exp 


/  a2  +  /?2  \ 
\  2W  ) 


iFi 


2  WiFi 


(4.59) 


where  T(2)  =  1  since  F(n)  =  (n  —  1)!.  From  Eq  (  4.40)  in  Section  4.1.3,  the  hypergeometric 
function  in  the  above  equation  reduces  to  the  form  (1  —  z).  Therefore,  the  expected  value  of  Qa+n 
is 


E(Q}+n) 


2W 


a2  +  {P  \ 

21 V  ) 
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E(Qs+„ )  =  2W  +  a2  +  /32 


(4.60) 


Now,  the  expected  value  of  Q^+n  is  derived  here.  From  Eq  (  4.53),  the  mean  square  value  is 


/OO 

Q^+nP{Qs+n)dQs+n 

*oo 

=  fir”-  (ffivs+y) 

■  i?esp('£L2 %r)l 


where 


i  =  y/Q^n 


dt  — 


l 


dQ*+» 


As  before,  from  Eq  (  4.31)  the  expected  value  of  Q%+n  is 


£(<??+„)  = 


r(3)  (wf 

— — —exp 


31 VVi 


a2  +  /?2\ 

2W  ) 


(4.62) 


The  confluent  hypergeometric  series  in  the  above  equation  can  be  expanded  to  a  finite  series  using 
Eq  (  4.32)  from  Section  4.1.3.  This  expansion  is  as  follows: 


,p.  (-2  i-s'j  =  i  _  j.  _ 

A  1  \  )  }  /  w  1  ^ 


(4.63) 


This  expansion  leads  to  the  following  form  of  the  mean  square  value  of  Q3+n: 
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E(QUn)  =  8W2 


Q2+/?2  (a2  +  /?2)2 

W  +  8iy2 

=  8  W2  +  8 W  (a2  +  P2)  +  (a2  +  /?2) 2 


(4.64) 


The  last  expectation  to  derive  is  the  autocorrelation  function  of  the  output  of  a  quadratic 
detector.  The  autocorrelation  function  is  given  by 

/oo  roo 

••■/  (x2  +  y2)  (xl  +  yl)p(xi,yl;x2,y2)dxidyldx2dy2  (4.65) 

■CO  j  —  CO 

where  p(^i,  j/i; 1/2)  is  given  by  Eq  (  4.46).  Lawson  and  Uhlenbeck  provide  the  solution  to  these 
integrals  (11:155).  The  autocorrelation  function  of  the  output  of  a  quadratic  detector  resolves  to 


Q1Q2  —  (cvj  +  $1  +  2TV)  {ot\  $2  +  2W) 

+4W>(c*i<*2  +/?i/?2)  +  4l72p2 


(4.66) 


V*.  Effects  of  CW  Interference 


In  this  chapter  the  effects  of  single  frequency  or  continuous  wave  (CW)  interference  are 
explored.  First,  the  effects  of  CW  interference  on  a  linear  detector  is  analyzed.  Finally,  the  effects 
on  a  quadratic  detector  are  determined.  This  analysis  closely  follows  the  analysis  presented  by 
Lawson  and  Uhlenbeck  (11:347-353).  For  this  analysis,  a  square  radar  pulse  is  assumed  with  pulse 
width  r  and  amplitude  5.  The  pulse  is  defined  by  the  expression 


s(t)  = 


Ssin  (2  nf,t  +  6,) 
0 


M<i 

elsewhere 


(5.1) 


where 

/*  =  signal  frequency 
6S  =  random  signal  phase 

It  is  further  assumed  that  the  signal  frequency,  /,,  is  centered  in  the  IF  filter  passband. 
The  CW  interference  is  represented  by  E(t): 


E(t)  =  £7  sin  (2nfcwt  +  0 

cu/) 


(5.2) 


where 

fcw  =  CW  interference  frequency 

9CW  =  random  ase  of  the  CW  interference 

It  is  assumed  that  the  voltage  of  the  CW  interference,  £7,  is  much  greater  than  both  the  signal 
voltage,  5,  and  the  average  value  of  the  noise  voltage.  The  input  to  the  IF  filter  (neglecting  thermal 
noise)  is  the  sum  oLthe  pulsed  signal  and  the  CW  interference.  Let  st'n(*)  be  the  IF  filter  input: 


5-1 


sin  (0  — 


(5.3) 


Ssm(2irftt  +  0t)  +  Esm(2nfcwt  +  Ocw)  ,  |<|  <  § 

< 

Es\n(2i:fcwt  +  ecw)  '  ,  elsewhere 

It  will  be  convenient  to  put  s,-„(t)  into  quadrature  form  in  terms  of  the  center  frequency,  fa,  so  that 

Sin(t)  =  a0{t)  cos  2  it  fit  +  p0{l)  sin  2nf,t  (5.4) 

For  |<|  <  §,  S{n(t)  can  be  manipulated  using  trigonometric  identities  into  the  quadrature  form: 

Si„(t )  =  5(sinw,i  cos  d3t  +  cosw3  sin  0, ) 

4 

+E  [sin  (w'  +  )  t  cos  0CW  +  cos  (w'  +  w, )  f  sin  0CU)]  (5.5) 

where 

=  27 xfs 

U)f  =  wCK  -  Wi 

Sin(t)  =  S(s\ncjstcos09  +cosiostsinOs) 

+E  [(sin  u /t  cos  ust  +  cos  u>*t  sin  u>3  J)  cos  0CW 
+  (cos  uft  cos  u>3J  ~sinu/Jsinu>3i)sin0Cu/] 

=  cosw3/[Ssin03  +  i?(sinu/*cos0cw  +  cosw^sin^)] 

+  sin  u;3i[S  cos  03  -f  E  (cosu'tcosOcw  -  smut* t  sin  0ctv)] 

=  cos  w3 1  [5  sin  0S  +  E  sin  (uft  +  Ocw)] 

+  sin ust  (5  cos  0S  +  Ecos(u'i  +  0eu,)]  (5.6) 
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Likewise,  when  \t\  >  -  the  input  signal  becomes 


Sin(t)  =  cosojstEsm(u)ft  +  0CW)  $\nu>stE cos(w*t  -t-  0CW) 
From  Eqs  (  5.6)  and  (  5.7)  a0(t)  and  f30(t)  can  be  expressed  as 


(5.7) 


m = { 


E  sin  [2;r  (fcw  -  /„ )  i  +  Qcw]  +  5  sin  0, 

.  |<l<§ 

E sin  [2<r  ( fcv .  -  f,)t  +  0iW\ 

,  elsewhere 

E  COS  [2 TT  (fcw  —  f»)t  +  0CU>]  +  S  cos  0, 

.  w<§ 

E cos  [2tt  (few  ~  ft)  ^  "i*  @cw\ 

,  elsewliere 

(5.8) 


(5.9) 


The  analysis  by  Lawson  and  Uhlcnueck  assumes  that  the  IF  filter  has  a  Gaussian  passband. 
The  transfer  function  of  this  filter  is  expressed  as 


H(f)  =  exp  (~%r)  (5-10) 

where  a  =  1.18,  so  that  the  half  power  bandwidth  of  the  filter  is  B.  Here,  the  frequency  will  be 
taken  about  the  frequency  /  =  which  is  the  center  frequency  of  the  IF  filter  passband.  The 
transfer  function  of  this  filter  is  illustrated  in  Figure  5.1  for  positive  frequencies. 

The  output  of  the  IF  filter,  sout(t)  can  be  calculated  by  Fourier  transform  techniques  as 
follows 


*>««(*)=  r (5.11) 

where  5>n(f)  is  the  Fourier  transform  the  input  to  the  IF  filter.  Sm(/)  can  be  manipulated 
through  properties  of  the  Fou  .  r  transform  t.  be 
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Figure  5.1.  IF  Filter  Transfer  Function 


Sin(f)  =  \  [A(f  -  fs)  +  A(f  +  /,)]  +  1  [B(f  -  / ,)  -  B(f  +  /,)]  (5.12) 


where 


(x0(i)e-i2*!tdt 

p0{t)e-i2*Sldl 


(5.13) 

(5.14) 


are  the  Fourier  transforms  of  a0  and  f30  respectively.  Using  Eq  (  5.12),  soui(i)  can  be  expressed  by 
the  following  equation: 
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(5.15) 


sout(t)  =  \J*  H(f)[A(f  -ft)  +  A(f  +  f,)]  efWdf 

/I  H(f)  [Bif  “  fs) "  B(/  +  /s)1 


This  expression  can  be  manipulated  to  the  form 


sout  CO 


/co  roo 

A(f)H(f)e^2lt^  df  +  sin  /  B(f)H(f)e^( df  (5.16) 

•CO  J  —  CO 

a(0cos27r/5t  -f  (3(f)  sin 2nfst  (5.17) 


where 


/oo 

A{f)H{fy2*l'df  (5.18) 

■CO 

/oo 

B(f)H(fy2*J(df  (5.19) 

■CO 


By  substituting  Eqs  (  5.8)  and  (  5.9)  into  the  above  equations  a(0  and  j3(t)  become  (11:348) 


a(t)  =  { 


m  =  < 


£?exp 

< 

Eexp 

«*(/cw-/.)_a  I 
B 3 

r 

Eex  p 

< 

E  exp 

S2(/ew7/»)3  j 

sin  {Jcw  -  fs)i  -t-  ucw\  -J-  o sin  vs  ,  |t|  <.  f 

sin  [2tt (few  -  fs)i  +  0CW]  ,  elsewhere 


,  elsewliere 


(5.20) 


(5.21) 
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5.1  Effects  of  CW  Interference  on  Linear  Detection 


In  this  section,  the  effects  of  single  frequency  CW  interference  on  linear  detection  are  analyzed. 
Since  the  CW  voltage  was  assumed  to  be  greater  than  either  the  signal  or  mean  noise  voltage,  the 
average  output  of  the  linear  detector,  Lt  when  the  input  is  sout(i)}  is  found  from  the  approximation 
of  Eq  (  4.37)  in  Section  4.1.3  to  be 


L  -  \fa2  +  /?2  (5.22) 

Using  the  expressions  for  a(i)  and  /?(i),  Eqs  (  5.20)  and  (  5.21),  the  average  value  for  \t\  <  § 
becomes 


L  =  [. X 2  sin2(a/t  4  0CW)  +  2X5 sin (w't  4*  0CW) sin 0S  4  52  sin2  0S 

4X2  cos2(o/J  4-  0CW)  4  2X5  cos(coft  4  0CW)  cos  0S  4-  52  cos2  09\ 1 
=  {X2  4-2X5  [sin(a/t  4  0CW)  sin  0S  cos  (u't  4  0CW)  cos  0S]  +  52}3  (5.23) 


where 


0)*  —  2/T {few  ~~  fa) 

X  =  E exp  |-23&g£&£] 


The  above  expression  for  L  can  be  reduced  through  trigonometric  identities  to 


L  =  [X2  +  2X5cos(w'<  +  <p)  +  S2]  ’ 


(5.24) 


where 


<P  =  0Cw-  0 , 
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It  was  assumed  that  the  amplitude  of  the  CW  interference,  E)  was  much  greater  than  the  amplitude 
of  the  signal,  S.  Therefore,  the  S 2  term  in  the  above  equation  is  negligible  and  L  can  just  as  easily 
be  written 


L  =  [. X 2  +  2  XS  cos  (a/*  -f  tp)  -f  S2  cos  2(u't  +  y?)]  ^ 

=  X  +  S  cos(a /t  -f-  <p) 


(5.25) 


When  |<|  >  S  =  0  and  L  =  X.  This  leads  to  the  following  expression  for  L: 


L  = 


j  2? exp  +  5 cos  [27r  (/cu,  -  /*)*  +  <?]  ,  |t|  <  \ 


Eex p 


,  elsewhere 


(5.26) 


Lawson  and  Uhlenbeck  assume  that  the  video  bandwidth,  6,  is  very  narrow,  so  that  only 
the  variation  inside  of  the  pulse  width,  r,  needs  to  be  considered  (11:349).  Therefore,  L  is  time 
averaged  over  r  to  yield  L.  For  the  case  of  the  signal  present,  5^0 


X  +  S  cos  (u't  +  ip)dt 


=  X  + 


oj’t 


.  j 

f  w'r  \ 

i  •  / 

sin  | 

[fT +”J 

|  -  sin  | 

|  +sin  | 

sm  | 

i-h 

u't  \\ 

--*)] 


(5.27) 


where 


w'  =  2t r(/CU)  -  ft) 
X  =  Eex  p 
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Without  the  signal  present,  5=0  and  the  above  equation  reduces  to  L  =  X.  With  this,  the 
equation  for  L  can  be  written 


Lcw+a  -  Eexp  J  L]  +  2^(/c«-/,)r 

L  =  i  x  {sin  [tt  (/cw  -  /4)  r  +  y>]  +  sin  [it  (fcw  -fs)r-<p]}  ,  with  signal  (5.28) 

^  Lcw  =  i?exp  j-  — 1 ^gf^- j  ,  without  signal 

To  determine  the  effects  of  the  CW  interference,  the  deflection  criterion  presented  by  Lawson 
and  Uhlenbeck  is  applied.  This  criterion  states  that  a  detection  has  occurred  if  the  average  value 
of  the  deflection  for  the  signal  plus  noise  over  the  average  value  of  the  deflection  for  noise  only  is 
on  the  order  of  magnitude  of  the  standard  deviation  of  the  deflection  for  noise  only  (11:161).  This 
criterion  in  equation  form  is  (11:163) 


h  :  =  k  (5.29) 

yy2N-(yN)2 

where  k  «  1.  For  the  case  of  the  linear  detector,  Ls+ at  is  derived  by  assuming  several  observations 
are  averaged,  and  that  the  phase,  ipy  is  uniformly  distributed  over  the  interval  (0,  2tt).  Further, 
a  central  limit  theorem  argument  allows  the  use  of  a  Gaussian  density  for  a  large  number  of 
observations  of  the  deflection.  The  mean  of  this  Gaussian  density  is  Lcw+S  —  Lcw  from  Eq  (  5.28) 
and  the  variance  is  2 W.  Since  Lcw+s  —  Lcw  will  vary  between  positive  and  negative  values  due  to 
phase  fluctuations,  the  average  will  be  taken  over  one-half  of  the  Gaussian  density.  As  Lawson  and 
Uhlenbeck  point  out,  Ls+n  is  defined  by  (11:349) 


2  f 2 *  1  /°° 

S+N  ”  V2rtV  Jo  2tt  J0  X 


exp  < 


~  1 2  \ 


®  {Lcw+s  ~  ^•'eu;)j 

w 


>  dxchp 


(5.30) 
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The  derivation  of  Ls+n  follows.  The  LCw+s  -  term  can  be  reduced  through  trigonometric 
identities  to 


&CW+9  ~  Lew  — 


”  /a)7* 
S 


{sin  [7 r  (fcw  -fs)r  +  (p]+  sin  [tt  (fcw  -  fs)r-<p]} 


w  —  fs)r 
=  7  cos  v? 


sin  [tt  (fcw  -  /A)r]cos<p 


(5.31) 


where 

*(fj-7;yrs[n{* (few  h)r ) 
Substituting  the  above  expression  into  Eq  (  5.30),  gives 


_  ^  r  2n  1 00 

Ls*n  =  7vml  L  *“p 


(a;  —  7  cos  v?)2 
ZW 


dxd<p 


1  f 


r  2tt  rOO  .  _  . 

- - :  /  /  ( uV2W  +  7  cos  tp )  e~u3V2Wdudip 

xVwvJo  > 

1  (  f°°  _  2  2  \ 

— 7=  /  [  /  ViWue^^du  4-  /  7cos^e“u2du  )  d(p 

Vy/7r  Jq  y  ^llv  ) 


(5.32) 


where 


.  _  X-7C0S  (p 

u-  vw 


du 


72 W 


The  component  integrals  of  the  above  equation  will  be  solved  here.  The  first  integral  reduces  to 


/OO  /  rOO  />  0 

y/2Wue~u2 du  =  y/2W  (  /  ne"tl2du  4-  / 

■  \Jo 

W  \  y/2  W 


ue~u  du 


5-9 


(5.33) 


The  second  integral  is  nearly  of  the  form  of  the  co-error  function,  and  can  be  written  in  terms  of 
the  co-error  function  by 


L 


7  cos  (pe~u  du  = 


COS( 


2 

tv'*7 


■  cos  <p  erfc(— z) 


(5.34) 


where 


*  —  7  cosy 

z  i/w 


This  expression  can  be  reduced  to  a  form  containing  the  error  function  by  the  use  error  function 
relationships  found  in  the  Handbook  of  Mathematical  Functions  as  follows  (13:297): 


erfc(-2)  =  1  —  erf(— 2) 
erf(-2)  =  -erf(z) 

erfc(— 2)  =  1  -f  erf(^) 


With  the  above  relationship,  the  integral  can  be  written 


7  cos  <pe  u  du  = 


7V^r 


2 


cos  tp  [1  4*  erf(2)] 


(5.35) 


(5.36) 
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A  series  expansion  for  erf(z),  from  the  Handbook  of  Mathematical  Functions  is  substituted  for  the 
error  function  above  (13:297): 


ycosipe~u  du 


2 


COS  (p 


2 


COS  tp 


(-l)nZ2"+l 
n!(2n  +  1) 

/  \ 2n+l 

(-1)"  (w) 

n!(2n  -f  1) 


(5.37) 


Substituting  Eqs  (  5.33)  and  (  5.37)  into  Eq  (  5.32),  Ls+n  can  be  written 


Ls+n  = 


.  iV* 

H — —  cos  (p 


.  ,  if  H>"  (^) 

V*  ^  n\(2n  +  1) 


2n+l 


nsO 

1  j  f2*  [W  (  72cos2^  j 

57?  U  VT  “f(— w-)1'*' 

[2*  7%/^  i  f27r  [7  cos  (p 

+  /o 


>  dip 


73  cos3 
G(2W)§ 


+  ■ 


(5.38) 


The  argument  of  the  first  integral  in  the  above  equation  can  be  expanded  in  a  power  series  of  the 
form  of  Eq  (1.211,3)  of  Gradshteyn  and  Ryzhik  (8:22).  With  this  series,  the  first  integral  becomes 


Keeping  only  the  terms  up  to  72,  the  integral  becomes 
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(5.40) 


The  second  integral  in  Eq  (  5.38)  reduces  to  zero.  By  keeping  only  the  terms  up  to  72,  the 
third  integral  becomes 


/*2  ir 

J  0 


7  cos  <p 


7  cos  <p  73  cos3  ip 


y/2W  6(2W)% 


d(p 


-L 


2  T  9  9 

Y  COSJ  Ip 


vW 


72tt 

vW 


(5.41) 


By  substituting  the  integral  results  of  Eqs  (  5.40)  and  (  5.41)  into  the  equation  for  Ls+n,  Eq 
(  5.3S),  Ls+n  becomes 


Ls+n  = 


1 

7T>/? 


/2W  *  i  .  il 

7T  |  W  L  x(fcw-f,)r 


sin7r(/ctu  ~/,)r 


(5.42) 


The  term  Tat  is  derived  by  setting  S  =  0  in  the  above  equation,  so  that 


Ln  = 


2W 


(5.43) 


Now,  Ls±n  —  Ln  is  obtained: 


-  -  s2  lm 


sin  n(fcw  -  fs)r 


^  ( few  ~  ft)  T 


(5.44) 
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Because  a  very  narrow  video  bandwidth,  6,  was  assumed,  the  above  equation  is  reduced  by  a 


factor  of  br  to  arrive  at  the  video  output: 


-  -  S2br  IW 


sin  n(fcw  -fs)r 

,  ^  ( few  “  fs)  T 


-  2 


(5.45) 


To  arrive  at  the  standard  deviation  of  the  CW  and  noise  only,  the  power  spectral  density  will 
first  be  obtained.  The  PSD  of  the  CW  plus  noise  is  derived  by  taking  the  Fourier  cosine  transform 
of  the  autocorrelation  of  the  output  of  the  linear  detector.  Since  it  was  assumed  at  the  outset  that 
the  CW  power  was  much  greater  than  the  noise  power,  the  large  signal  approximation  from  Section 
4.1.3,  Eq  (  4.49),  is  used  for  the  autocorrelation  function.  As  Lawson  and  Uhlenbeck  point  out, 
“since  E  ^  W,  the  main  contribution  to  the  continuous  spectrum  of  the  linear  detector  output 
come  from  the  term  proportional  to  p  ”  (11:350).  So,  the  autocorrelation  becomes 


LiL2  =  Wp(r) 


[(a?  +  /?12)(‘>i  +  $)}* 


(5.46) 


where  c*i  and  a  2  are  from  Eq  (  5.20)  with  5  =  0  at  times  <1  and  <2  respectively.  Likewise,  the  /? 
terms  are  from  Eq  (  5.21)  with  S  =  0.  Now,  the  autocorrelation  becomes 


r  r  ■  .  .  X2sin©i  sin  ©2  +  X2  cos©i  cos  ©2 

LiL2  =  Wp{T) - r 

[{X2  sin2  ©1  +  X 2  cos2  ©1)  (X2  sin2  ©2  +  X 2  cos2  ©2)] 3 
=  Wp(r)  (sin©i  sin  ©2  +  cos©i  cos  ©2) 

=  W p{r)  cos(©i  -  ©2)  (5-47) 
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where 


X  =  2?  exp  [~2!!£gfZi£J 

©1  —  2ir {few  ~  fs)i i  "4*  0Cw 
©2  =  2?T  (/cu,  -  /,)  i2  +  Ocw 

With  the  above  relationships,  the  autocorrelation  becomes 


LiL2  =  Wp{r)  cos  [2?r  (/cu,  -  /,)  U  -  2-rr  ( fcw  -  f,)t2] 


=  W p(f)  cos  [2?r  {few  -  f, )  r) 


(5.48) 


Next,  the  normalized  correlation  coefficient,  p(r),  is  calculated  by  the  use  of  Eq  (  4.45). 


(  s  fZo\H(f)\2cOS2*fTdf 
m  ■  fZ  \H{f)M 


(5.49) 


The  IF  filter  was  assumed  to  have  a  Gaussian  passband  whose  transfer  function,  H(f)}  was  defined 
by  Eq  (  5.10).  With  this  equation,  the  numerator  in  the  equation  for  p(r)  is 


J  \H (/)|2  cos 2tt frdf  =  2  J  exp  cos  27 rfrdf 


(5.50) 


This  integral  is  of  the  form  of  Eq  (3.896,4)  from  Gradshteyn  and  Ryzhik,  and  resolves  to  (8:480) 


r 


i//(/)  j2  cos  2tt / rdf 


„  fi  R& 

2\2V2?eXP 


(2<r  r)2  jB2  1 
4  (2a2)  Jj 


B  1 

l  i 

(  tttB  \  2 

O  eXP  “o 

a  V 

2  2 

\  a  J  _ 

(5.51) 
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The  denominator  reduces  to 


/: 


\m?df  = 


-  Kf)(^) 

B  fir 
a  V  2 


(5.52) 


Combining  Eqs  (  5.51)  and  (  5.52),  the  normalized  correlation  coefficient  becomes 


p(r)  =  exp 


1  /  7T  t\ 

~2  V“a 


:rrl?\  2 


With  the  above  the  autocorrelation  function  becomes 


(5.53) 


L1L2  =  Rr,{T)  —  W  exP 


1  / ktB V 


cosw  r 


(5.54) 


where 


U1  =  2x  (fcw  -  fs) 


The  power  spectral  density  is  calculated  by  taking  the  Fourier  cosine  transform  of  the  auto¬ 
correlation  function: 


SiU) 


A00 

2  /  /f^(r)coso;rdr 
Jo 


cos  a/r  cos  urdr 
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w  r  [ 1  f*rBY]  <  <  ^  a 

W  J  exp  I— -  J  cos(w  -wjrar 
+W  J  exp  — ^  (~a~)  cos(w'  +u)rdr 


The  integrals  in  the  above  equation  are  of  the  form  of  Eq  (3.896,4)  found  in  Gradshteyn  and  Ryzhik 
(8:480).  Using  this  equation,  the  PSD  reduces  to 


W 

2  V  ?r252 


f  r  aV-w)*l  r  aV  +  w)2!) 

Lexp  [— W^-J  +  [ — } 


Wa  /Tf 

"  25  V  ir  | 


[  2a2  (/cm  —  /»  +  /)2 


This  PSD  is  illustrated  in  Figure  5.2. 


Since  a  narrow  video  bandwidth,  6,  was  assumed,  the  variance  of  the  CW  plus  noise  video  is 
2Sx(0)6.  Therefore,  from  Eq  (  5.56)  the  variance  of  the  CW  plus  noise  video  output  is 


nLWa  1 

T  f 

2blB\l 

*{eXp 

2  aWb  1 

T 

B  V 

—  exp  ■ 

7T 

2a2  (fcw-f.? 


B2 

2a2  (/cu,-/,)2 

B2 


+  exp 


2a2  (/cu,-/,)2 

B2 


(5.57) 


The  standard  deviation  of  the  video  output  is  simply  the  square  root  of  the  variance.  So,  the  square 
root  of  Eq  (  5.57)  is 


aHfcW-f,)2 


B2 


(5.58) 


From  the  equation  for  the  deflection  criterion,  Eq  (  5.29),  an  expression  for  the  minimum 
detectable  average  signal  powe*  is  derived.  The  deflection  criterion  is  represented  by 


lI-{ln)2}* 


(5.59) 


where  C  is  a  constant  determined  by  the  radar  receiver  structure.  Substituting  Eqs  (  5.45)  and 
(  5.58)  into  the  above  equation  yields 


4  W 


w 

sin7r(/CU)  —  f,)r 

2  -c(2)\l2aWbexo 

a2  (few  ~  f>)2 

V  7T 

.  * (fcw-fs)r  . 

v  b  exp 

~B 2 

(5.60) 


Lawson  and  TJhlenbeck  define  the  minimum  detectable  average  signal  power  as  Bm\n  = 
where  T0  is  the  pulse  recurrence  interval  (11:350).  Now  /Jmjn  is  derived: 
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min 


S2t 


To 

4  CW 
bT0 


fabir  ( 2\* 

'  a2  {few  —  ft)2 

’  ft  (few  —  fs)  T  ' 

V  B  \x)  6XP 

B2 

s\mr(fCw  -f»)r. 

(5.61) 


From  Eq  (  4.4)  of  Section  4.1.1,  W  can  be  written  in  terms  of  the  IF  bandwidth,  B,  by 


/oo 

=  <W> 


Substituting  the  above  result  into  Eq  (  5.61)  the  minimum  detectable  average  signal  power  becomes 


min 


4 <t2CB  HF 
abT0  \J  2 

4  a2C 


labn  f  2\  * 

'  a2  (few  —  /s)2 

ft  (few  —  fa)  ft  ' 

V  B  [*J  p 

B2 

sinir (few- f»)T. 

inn  . 

VS5(  2')'"p 


«2  {few  ~  ft)2 


B2 


ftjfew  -  f»)r 


sin7r (fcw  -fs)r J 


(5.63) 


To  quantify  the  effects  of  on-frequency  CW  interference,  the  minimum  average  power  required 
for  the  case  of  large  detuning,  Pm  in(°°)>  is  first  derived.  Large  detuning  implies  that  |/cty-/3  |  >  B , 
and  in  this  case  the  CW  interference  can  be  considered  to  be  nonexistent  (11:351).  This  derivation 
°f  ^min(°°)  rnirrors  that  of  the  derivation  for  the  minimum  signal  power  when*  the  GW  was  present. 
First,  the  numerator  of  the  expression  for  the  deflection  criterion,  Ls+n  —  Ls+n  is  derived.  Ls+n 
is  obtained  from  the  equation  for  the  expectation  of  Eq  (  4.34)  in  Section  4.1.3: 
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Next,  the  denominator  for  the  deflection  criterion,  ^L2S  —  (LN)2j  ^  is  obtained  from  the  PSD 
for  the  noise  only  case.  First,  the  Fourier  cosine  transform  of  the  autocorrelation  function,  Rl(t) 
is  taken  to  obtain  the  PSD,  5x,(/).  Eq  (  4.48)  gives  the  small  signal-to-noise  ratio  approximation 
for  L1L2 .  For  the  noise  only  case,  5  =  0  and  a(t)  =  fi(t)  =  0.  Therefore,  the  only  applicable  term 
of  this  approximation  is  the  term  proportional  to  p2,  and  the  autocorrelation  function  becomes 

ThLi  =  Rl{t)  =  *Wp*ST\  (5.70) 

The  normalized  correlation  coefficient,  p(r),  was  defined  by  Eq  (  5.53).  Taking  the  Fourier  cosine 
transform  of  the  above  t  'nation  yields  the  PSD: 


/•OO 

Si(/)  =  2  /  Rl(t)  coswrdr 

_  ,  r*w?<T) 

Jo  8 

4*  r  f  /  jtt. 

=  wL  exp  \-{- 


cos  'jjrdr 


cos  lot  dr 


(5.71) 


The  above  integral  is  of  the  form  of  Eq  (679)  from  the  CRC  Tables  and  reduces  to  (16:466) 


SlU)  = 


Wab^/K 


8  B 


*  (  a2/2\ 


(5.72) 


As  before  for  a  narrow  video  bandwidth,  the  variance  is  25z,(0)6: 


Ll-(LNf  = 


4  B 


(5.73) 


Substituting  Eqs  (  5.69)  and  (  5.73)  into  the  expression  for  the  deflection  criterion  one  finds 
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S2br  fnW  „(  Waby/i r  * 
AW  V  2  ”  \  AB  ) 


(5.74) 


where  W  was  defined  by  Eq  (  5.62).  This  equation  can  be  reduced  to  the  expression  for  ^min(°°) 
as  was  doV\e  for  the  case  with  CW  present: 


pmin(°°)  = 


A<t2C 

To 


(5.75) 


For  on-frequency  CW,  ( fcw  =  /,),  the  required  minimum  average  signal  power  relative  to  the 
power  required  when  no  CW  is  present  is  given  by 


jPmin(O)  = 


(8*2)* 


=  4.75  dB 


(5.76) 


Lawson  and  Uhlenbeck  compare  this  value  with  a  value  of  3.5  dB  obtained  experimentally  by  A. 
L.  Gardener  and  C.  M.  Allred  (11:351).  Lawson  and  Uhlenbeck  also  compare  experimental  data 
with  theoretical  calculations  for  various  values  of  r\few  —  f$\*  The  experimental  data  was  collected 
with  a  narrow  video  bandwidth  and  a  large  IF  bandwidth  where  Br  =  2.2  (11:351).  This  data  is 
displayed  in  Table  5.1.  As  Lawson  and  Uhlenbeck  point  out,  the  theoretical  values  of  infinity  for 
integer  values  of  r\fcw  -  fs\  arise  because  S' 2  terms  were  neglected  (11:352).  Further,  theoretical 
values  for  r\fcw  -  fs  \  >  2.5  were  not  calculated  since  the  equation  for  Pm\n  is  not  valid  for  values 
of  | few  —  fs\>  B  (11:352).  The  theoretical  values  in  Table  5.1  are  calculated  from  the  following 
equation: 


P  . 

1  min 

-^min(^) 


=  exp 


a2(fcw~f,)2 

( few  ~~  fs)  T 

B7 

sin  v{few-f»)r. 

(5.77) 
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Table  5.1.  Theoretical  vs  Experimental  Values  of  Signal  Threshold  Due  to  CW  Detuning  for  a 
Linear  Detector 


Akw-fA 

PrmJpmm(0)>  dB 

Experimental 

Theoretical 

s  0.5 

5.5 

3.6 

i  1 

13.5 

00 

1.5 

11 

10.75 

2 

15 

00 

2.5 

9 

10.3 

3 

5 

3.5 

0 

1  4 

-2.5 

4.5 

-3.5 

(11:352) 


5.2  Effects  of  CW  Interference  on  Quadratic  Detection 

In  this  section  an  analysis  of  the  effects  of  CW  interference  on  quadratic  detection  is  accom¬ 
plished.  The  analysis  here  closely  follows  that  of  the  analysis  of  the  effects  on  linear  detection  of 
the  previous  section.  The  average  output  of  the  quadratic  detector  for  the  case  of  signal  plus  noise 
is  given  by  Eq  (  4.60)  of  Section  4.2.1: 


Q  =  2W  +  a2  +  02  (5.78) 

Substituting  the  expressions  for  a(t)  and  /?(*)  ,  Eqs  (  5.20)  and  (  5.21),  the  average  value  of  the 
detector  output  for  \t\  <  |  becomes 


Q  =  2 W  +  X2  +  6CW)  +  2XS s\n(ojft  +  Ocw)sinOs  +  S2sin20, 

+X2  cos  2(w,t  +  9CW)  +  2XS  cos  (a H  -f  0OU,)  cos  6S  +  S2  cos2  9S 
=  2 W  +  yY2  +  2 XS  [sin (w't  +  0CW) sin 9S  cos (uft  +  0CW)  cos  0S\  +  S2  (5.79) 
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where 


u*  —  2n(fcw  —  /,) 

X  =  Eex p  [-aa(^r/,)2] 

By  the  use  of  trigonometric  identities,  the  above  expression  for  Q  can  be  reduced  to  the  following 
form: 

Q  =  2W  +  X2  +  2  XS  cos(u H  +  <p)  +  S 2  (5.80) 


where 

<P  —  0CW  —  0S 

Since  it  was  assumed  that  E  S,  the  S2  term  in  the  above  equation  is  negligible.  Therefore, 
neglecting  the  S2  term,  Q  can  be  written 


Q  =  2W  +  X7  +  2XS  cos(u H  +  *>)  (5.81) 

When  \t\  >  S  =  0  and  Q  =  2W  +  X 2.  This  leads  to  the  following  expression  for  Q: 


Q-\ 


2W  +  E2  exp  [-2<fl/gs=/«E] 

+2i?S  exp  f— a^cgr^*^]  cos  [27r  (/cu,  -  /,)  t  +  v>] 
2W  +  E2  exp  [-2o-(/^,-/,)-] 


i*i<§ 

elsewhere 


(5.82) 


For  a  very  narrow  video  bandwidth,  b ,  only  the  variation  inside  of  the  pulse  width,  r,  needs 
to  be  considered.  Therefore,  Q  is  time  averaged  over  the  pulse  width,  r,  to  arrive  at  Q.  For  the 
case  of  the  signal  present,  5^0 
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1  fa 

=  -  /  2W  +  X2  +  2XScos{w't  +  <p)dt 

T  J_  J 

v,  2XS  r .  (w't  \  .  i 

f  w't  \1 

=  2W  +  X2  +  —  [sm^—  +  ¥>J-sin| 

r-+v)J 

2 XS  f .  (w't  \  .  , 

rwv  m 

=  2W  +  X2+  u,r  |sm  ^  2  +<p)+  sin  | 

^  2  *7] 

(5.83) 


where 

u'  =  2n(fcw  -  /,) 

X  =  Eex p 

When  no  signal  is  present,  5  =  0,  and  the  above  equation  reduces  to  Q  =  2 W  +  X 2.  With  these 
expressions,  the  equation  for  Q  becomes 


Qc*+)  =  21 V  +  E2  exp  p »-%-/■«£]  +  7(ycffyt>  exp 
Q  =  ^  x  {sin  [7T (/cu,  -  /i)  r  +  9]  +  sin  [ff  (/cw  -  /#)  r  -  <p]}  ,  with  signal  (5.84) 

Qcw  =  2W  +  i?2  exp  ,  without  signal 

As  was  the  case  for  the  analysis  of  the  linear  detector,  the  deflection  criterion  will  be  applied. 
As  before,  Qs+n  is  defined  by  (11:349): 


Qs+n  — 


2  /*2tt  ^  /*oo 


III 

Jo  to  Jo 


a;  exp  < 


(Qcw+j  Qcw)] 

2 W 


2  N 


>  dxdtp  (5.85) 


VtoW, 

Now,  Q^xiv  is  derived.  The  Qcu,x5.—  Qcu,  term  reduces  by  way  of  trigonometric  identities  to 
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Qew+t  Qcw  ~ 


SE 


n(fcw-ft)r 


exp 


«2(/cu)  ~  ft)* 


B2 


x  {sin  [tt  ( fcw  -  f, )  t  +  <p]  +  sin  [v  (fcw  -f,)r-<p]} 
2SE 


*( few  ~  ft)T 
—  7  cos  <p 


exp 


a2  {few  -  f*)2]  .  r  ft  t\  1 
— - 1  sin  [tt  {few  -  ft )  r]  cos  ip 


B2 


(5.86) 


where 

T  =  rJ/H-TTF  exP  [~  sin  I* (/c «  ~ft)r ] 

Substituting  the  above  expression  for  QCUI+,  -  Qcw  into  Eq  (  5.85),  gives  the  following  equation  for 
Qs+n: 


—  wbff- 


( x  -  7  cos  v?)2 


2  W 


dxdtp 


(5.87) 


The  integrals  in  the  above  equation  were  solved  in  Section  5.1.  They  reduce  in  terms  of  7  to  the 
form  of  Eq  (  5.42): 


Q*+*  =  \/v(1+4 Vi) 


(5.88) 


Substituting  the  expression  for  7  into  the  above  equation  yields  the  following  equation  for  Qs+n: 


■i/?{ 


S2E2 


Qs+n  =  \I  —  \i  +—<*p 


2a2(/cu)  ~/«)2l  [sin  7T  {fcw  -  ff)T 


B2 


[sin 

L 


1 2^ 


( few  /#)  ^ 


(5.89) 


The  term  is  derived  by  setting  S’  =  0  in  the  above  equation,  so  that 
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(5.90) 


Now,  Qs+n  —  Qn  is  obtained: 


IW  r 
VTexp[- 


2a2(/cu '-I.? 


B* 


sin  n(fcw  -  f,)r 


J  L  ^  ( few  f$)  T 


(5.91) 


Since  6,  the  video  bandwidth,  was  assumed  to  be  very  narrow,  the  above  equation  is  reduced 
by  a  factor  of  6r.  Therefore,  Qs+n  -  Qn  becomes 


Qs+n  ~  Qn  - 


S2EHt  [2W 

w 


exp 


2g2(/cu,  -  f,)2]  sin 7r(/eu,  -  f,)r 

B 2  J  [  *{fcw-f,)r  J 


12 


(5.92) 


To  arrive  at  the  standard  deviation  of  the  CW  and  noise  only,  the  power  spectral  density 
will  first  be  obtained  from  the  Fourier  cosine  transform  of  the  autocorrelation  function.  The  auto 
correlation  function  is  given  by  Eq  (  4.66)  of  Section  4.2.1.  Since  the  only  term  contributing  to 
the  continuous  spectrum  of  the  output  is  that  term  proportional  to  p,  the  autocorrelation  becomes 


Q1Q2  =  4 Wp(r)  (a\ot2  +  P1P2)  (5.93) 

where  c*i  and  a 2  are  from  Eq  (  5.20)  with  5  =  0  at  times  ti  and  1 2  respectively.  Likewise,  the  /3 
terms  are  from  Eq  (  5.21)  with  5  =  0.  Now,  the  autocorrelation  becomes 


Q  Q2  =  4 Wp{r)X2  (sin  ©1  sin  ©2  +  cos  ©1  cos  ©2) 
=  AW p(r)X2  cos(©i  —  ©2) 


(5.94) 
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where 


X  =  E  expj— — 

©l  =  2?r  (/ct  —  ft)ti  + 

©2  =  2?r  (/CU)  —  ft)  t2  4*  $eu; 

With  the  above  relationships,  the  autocorrelation  becomes 


Q1Q2  =  AWE?  exp 
=  AWE2  exp 


2g2(/cu>  ~/,)2 

j B2 

2o2U*  -f*)2 

B- 


p(t)  cos  [2jr  (few  -  /,)<i  -  2tt  (/Clu  -  /,)  <2] 


p(r )  cos  [2tt  (/eu,  —  /»)  r] 


(5.95) 


The  normalized  correlation  coefficient,  p(r),  was  calculated  in  Section  5.1  by  the  use  of  Eq 
(  4.45): 


cW-«p[-l  (^) 


(5.96) 


With  the  above  equation,  the  autocorrelation  function  becomes 


Q1Q2  =  Rq(t)  =  AWE2  exp 


2a2{fcw- fs)2' 

[  1  fir  tB\2' 

B2  J 

CXP[  2V  a  ) 

cosw'r 


(5.97) 


where 


u'  =  27 r(/cu,  -/,) 


The  power  spectral  density  is  calculated  by  taking  the  Fourier  cosine  transform  of  the  auto¬ 
correlation  function: 
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JfOO 

f  Rq(t)cosut(It 
0 


Figure  5.3.  PSD  of  the  Output  of  a  Quadratic  Detector 


Because  a  very  narrow  video  bandwidth,  6,  was  assumed,  the  variance  of  the  CW  plus  noise 
video  output  is  2S,<?(0)6.  Therefore,  (from  Eq  (  5.99)  the  variance  of  the  CW  plus  noise  video  output 
is 


Ql-(QJ 


2  _ 


8  aWbE2 
B 


4 


exp  I- 


4a2  (few  —  f,)2 


B 2 


(5.100) 


The  standard  deviation  of  the  video  output  is  simply  the  square  root  of  the  variance.  So,  the  square 
root  of  Eq  (  5.100)  is 


[%-«„)’]  *  =  (!)' W^ir'xp 


I,  2a2  (fcw  -  f,)~ 
B2 


(5.101) 


From  the  equation  for  the  deflection  criterion,  Eq  (  5.29),  an  expression  for  the  minimum 
detectable  average  signal  power  is  derived.  The  deflection  criterion  is  represented  by 


N  .  =  C  (5.102) 

where  C  is  a  constant  determined  by  the  radar  receiver  structure.  Substituting  Eqs  (  5.92)  and 
(  5.101)  into  the  above  equation  yields 


S-E2bT 

w 


h  w 

2  a2(fcw-ftf 

’sin?r (/CU)  -  f,)r 

VTexp 

B2 

.  * (fcw-fs)r  . 

c 


(j;Y24 


2aWb 

B 


■exp 


2a2  {few  -  /»)2 


B2 


(5.103) 


The  minimum  detectable  average  signal  power  was  defined  as  where  T0  is  the 

pulse  recurrence  interval  (11:350).  Now,  Pm,n  for  the  quadratic  detector  is  derived: 


5-29 


mm 


52r 

r0 

2CW 

£6T0 


/a&jr 


ir 


(fcw-fs)r 


[sin7r(/cw  -  /j)rJ 


(5.104) 


W  was  derived  in  terms  of  the  IF  bandwidth  B  in  section  5.1,  and  is  given  by  Eq  (  5.62) 

where 


W  = 


<r2B  [if 

~V  2 


(5.105) 


Substituting  the  above  result  into  Eq  (  5.104)  the  minimum  detectable  average  signal  power  becomes 


mm 


2  ^CB 
EabT0 

2a2C  firB 
T0EV  2 ab 


&/¥(!)* 

(2»)* 


ft  ( few  ~  fs)  T  ] 


sin  7r(/cu,  -  f9)r  J 

7T(/cti;  ~/j)t  12 


sin  ft(fcw  -  fs)r 


(5.106) 


An  analysis  of  the  effects  of  the  CW  interference  similar  to  that  for  the  case  of  a  linear  detector 
is  accomplished  here.  ^mjn(oo)  is  derived  to  quantify  the  effects  of  on-frequency  jamming.  The 
numerator  of  the  expression  for  the  deflection  criterion,  Qs+N  —  QN  is  derived  first.  is 

obtained  from  Eq  (  4.60)  of  Section  4.2.1: 


Qs+„  =2W  +  a2  +  02  (5.107) 

From  the  expressions  for  a(t)  and  /?(*),  Eqs  (  5.20)  and  (  5.21),  a2  +  ft2  =  S 2  when  no  CW 
interference  is  present.  This  leads  to 
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Qs+n  =  2W  +  S> 


(5.108) 


The  term  QN  follows  directly  from  Eq  (  4.55): 

Qn  =  2W  (5.109) 

Combining  the  two  equations  above  yields  the  following  expression: 

Qs+n~Qn=S2  (5.110) 

The  above  equation  is  reduced  by  a  factor  offer  to  account  for  the  narrow  video  bandwidth.  This 
results  in 

Qs+N-^s=s2br  (5-111) 

Without  CW  interference  present,  the  denominator  of  the  expression  of  the  deflection  cri- 
terion,  [of,  _  (qn  )2j  ^  is  derived  from  the  PSD  for  the  case  of  noise  only,  Sq( /).  First,  the 
autocorrelation  function  is  obtained  by  the  application  of  the  expression  for  Q1Q2,  Eq  (  4.66)  of 
Section  4.2.1.  In  the  noise  only  case  5  =  0  and  a  =  /?  =  0.  Therefore 


QxQ2  =  Rq(t)  =  4  wy  (r)  (5.112) 

The  PSD  is  derived  by  taking  the  Fourier  cosine  transform  of  the  autocorrelation  function: 
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SqU)  = 


2  /  Rq{t)  cos uftdr 

Jo 

=  8W2  J  p2(r)coswrdr 

■  — r-f-w 


coswrdr 


(5.113) 


The  integral  in  the  above  equation  was  solved  in  the  previous  section  by  the  use  of  Eq  (679)  from 
the  CRC  Tables  (16:466).  It  reduces  to 


4‘. /2a  /  a2/2\ 


The  variance,  Q2N  —  {QN)2,  is  merely  2Sq(0)6: 


(5.114) 


8W2a& 


Q2  -(Q  )2  = 


(5.115) 


Substituting  Eqs  (  5.111)  and  (  5.115)  into  the  expression  for  the  deflection  criterion,  Eq 
(  5.29),  one  obtains 


<-> 

Similar  to  the  derivation  in  the  previous  section,  the  above  equation  reduces  to  the  following 
equation  for  Pmm(°°): 


(5.117) 


To  compare  on-frequency  CW  jamming  with  the  case  of  large  detuning  one  obtains 
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jPmin(O)  =  i/fV 

*W«>)  A2/ 

=  i(1.49)  (5.118) 

This  value  for  the  on-frequency  jamming  is  contrasted  with  the  value  of  4.75  dB  obtained  for  linear 
detection.  Here,  the  magnitude  of  the  CW  jamming  signal,  E)  has  a  direct  impact  on  the  minimum 
power  required. 

Although  no  experimental  data  was  available  for  the  quadratic  detector,  theoretical  values  of 
Pmin/^min(®)  are  derived  for  values  of  r\fcw  -  /j|  <  2,5.  These  values  are  tabulated  in  Table 
5.2.  As  before,  Br  —  2.2.  For  values  of  r\fcw  -  /,|  >  2.5,  the  large  detuning  causes  the  effects 
of  the  CW  interference  to  be  negligible.  This  follows  from  the  experimental  data  obtained  for  the 
linear  detector.  The  theoretical  values  in  Table  5.2  are  calculated  from 


^min  _  n(fcw  -  jj) 7 
?min(°)  Ls*nir(/ew  —  f»)rm 


(5.119) 


Table  5.2.  Theoretical  Values  of  Signal  Threshold  Due  to  CW  Detuning  for  a  Quadratic  Detector 


r\fcw-  f,\ 

^min/^min^)’ 

0.5 

3.9 

1 

oo 

1.5 

13.5 

2 

oo 

2.5 

17.9 
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VL  Conclusions  and  Recommendations 


6.1  Conclusions 

From  the  expression  for  the  minimum  power  required  for  signal  detection  by  the  linear  de¬ 
tector,  Eq  (  5.63),  the  minimum  power  is  a  function  of  the  frequency  difference  fcw  -  /,.  For 
on-frequency  interference,  fcw  —  Lawson  and  Uhlenbeck  show  from  experiment  that  the  mini¬ 
mum  required  signal  power  is  3.5  dB  higher  than  that  required  when  no  CW  interference  is  present 
(11:351-352).  This  increase  in  required  power  is  compared  with  the  theoretical  value  of  4.75  dB 
(11:352).  The  theoretical  and  experimental  data  in  Table  5.1  show  that  the  minimum  power  re¬ 
quired  for  detection  increases  as  the  detuning  increases  until  the  detuning  is  so  great  as  to  render 
the  CW  interference  negligible  (11:352). 

Similarly,  the  minimum  required  power  for  the  quadratic  detector  is  given  by  Eq  (  5.106). 
This  equation  shows  that  the  minimum  power  for  the  quadratic  detector  is  a  function  of  both 
the  frequency  difference,  fcw  —  /,,  and  the  magnitude  of  the  CW  signal,  E .  The  minimum  power 
required  when  the  interference  is  on-frequency  versus  the  power  required  without  CW  interference  is 
^(1.49).  The  data  tabulated  in  Table  5.2  show  that  as  is  the  case  for  linear  detection,  the  minimum 
power  required  increases  until  the  detuning  is  so  great  that  the  effects  of  the  CW  interference  are 
negligible. 

The  effect  of  a  pure  CW  tone  as  an  interfering  signal  is  not  a  significant  degradation.  The 
effect  of  CW  interference  is  especially  limited  when  the  frequency  of  the  CW  signal  differs  greatly 
from  that  of  the  carrier  frequency  of  the  pulsed  radar  signal.  However,  when  some  other  modulation 
is  present  with  the  CW  interference,  the  combined  effects  should  be  much  greater  than  that  of  an 
unmodulated  CW  signal. 


6-1 


6.2  Recommendations 

It  is  recommended  that  follow-on  work  be  accomplished  on  this  study  of  the  effects  of  noise 
jamming  on  radar  receivers.  This  follow-on  work  would  include  applying  the  analysis  techniques 
presented  here  to  the  effects  of  CW  interference  on  different  pulsed  radar  receiver  types.  Specifically, 
these  receivers  would  include  matched  filter,  pulse  compression  and  CFAR  receivers.  Additionally, 
various  jamming  waveforms  should  be  substituted  for  the  CW  interference,  and  the  effects  on  the 
various  receivers  should  be  analyzed.  Lastly,  instead  of  the  application  of  the  deflection  criterion, 
an  analysis  by  the  use  of  a  likelihood  ratio  test  should  be  used  to  determine  the  effects  of  the 
different  jamming  waveforms  on  the  radar  receiver’s  Pd  and  P/a. 


6-2 


Bibliography 


1.  Bennighof,  R.  H.  and  others.  “Effectiveness  of  Jamming  Signals.”  In  Boyd,  J.  A.  and  others, 
editors,  Electronic  Countermeasures,  Los  Altos  CA:  Peninsula  Publishing,  1978. 

2.  Cassara,  F.  and  others.  “A  Uniform  Power  Spectral  Density  Jamming  Signal,”  Proceedings  of 
the  IEEE ,  £7:331-332  (February  1979). 

3.  Davenport,  Jr.,  Wilbur  B.  Probability  and  Random  Processes ,  An  Introduction  for  Applied 
Scientists  and  Engineers.  New  York:  McGraw-Hill  Book  Company,  1970. 

4.  Davenport,  Jr.,  Wilbur  B.  and  William  L.  Root.  An  Introduction  to  the  Theory  of  Random 
Signals  and  Noise .  New  York:  McGraw-Hill  Book  Company,  1958. 

5.  DiFranco,  J.  V.  and  W.  L.  Rubin.  Radar  Detection .  Dedham  MA:  Artech  House,  Inc.,  1980. 

6.  Fitts,  LtCol  Richard  E.,  editor.  The  Strategy  of  Electromagnetic  Conflict .  Los  Altos  CA: 
Peninsula  Publishing,  1980. 

7.  Golden,  Jr.,  August.  Radar  Electronic  Countermeasures.  Washington  DC:  American  Institute 
of  Aeronautics  and  Astronautics,  Inc.,  1987. 

8.  Gradshteyn,  I.  S.  and  I.  M.  Ryzhik.  Table  of  Integrals,  Series,  and  Products  (Fourth  Edition). 
New  York:  Academic  Press,  Inc.,  1965. 

9.  Jordan,  J.  P.  and  P.  L.  Penney,  editors.  Transmission  Systems  for  Communications  (Fifth 
Edition).  Holmdel  NJ:  Bell  Telephone  Laboratories,  1982. 

10.  Knorr,  Professor  Jeffrey  B.  and  Maj.  DimitriosT.  Karantanas.  “Simulation  Optimizes  Noise 
Jammer  Design,”  Microwave  Journal ,  28: 273-277  (May  1985). 

11.  Lawson,  James  L.  and  George  E.  Uhlenbeck,  editors.  Threshold  Signals.  New  York:  McGraw- 
Hill  Book  Company,  1950. 

12.  Maksimov,  M.  V.  and  others.  Radar  Anti- Jamming  Techniques .  Dedham  M A:  Artech  House 
Books,  1979. 

13.  National  Bureau  of  Standards,  U.  S.  Department  of  Commerce.  Handbook  of  Mathematical 
Functions.  AMS  55.  Washington  DC:  Government  Printing  Office,  1964. 

14.  Papoulis,  Athanasios.  Probability,  Random  Variables,  and  Stochastic  Processes.  New  York: 
McGraw-Hill  Book  Company,  1965. 

15.  Schlesinger,  Robert  J.  Principles  of  Electronic  Warfare .  Englewood  Cliffs  NJ:  Prentice-Hall, 
Inc.,  1961. 

16.  Selby,  Samuel  M.,  editor.  Standard  Mathematical  Tables  (Twenty-first  Edition).  Cleveland: 
The  Chemical  Rubber  Co.,  1973. 

17.  Shanmugan,  K.  Sam  and  A.  M.  Breipohl.  Random  Signals  Detection ,  Estimation  and  Data 
Analysis.  New  York:  John  Wiley  and  Sons,  1988. 

18.  Skolnik,  Merrill  I.  Introduction  to  Radar  Systems  (Second  Edition).  New  York:  McGraw-Hill 
Book  Company,  1980. 

19.  Turner,  F.  M.  and  others.  “Noise  Quality  Optimizes  Jammer  Performance,”  Electronic  War¬ 
fare/Defense  Electronics ,  0:117-122  (November/December  1977). 

20.  Van  Brunt,  LeRoy  B.  Applied  ECM ,  Volume  1.  Dunn  Loring  VA:  EW  Engineering,  Inc.,  1978. 

21.  van  der  Ziel,  Aldert.  Noise.  New  York:  Prentice-Hall,  1954. 


BIB-1 


22.  van  der  Ziel,  Aldert.  Noise  in  Solid  Stale  Devices  and  Circuits.  New  York:  John  Wiley  and 
Sons,  1986. 


BIB-2 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hou'  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  Information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202*4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704*0188),  Washington.  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

December  1990  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE 

A  Study  of  the  Effects  of  Noise  Jamming  on  Radar  Receivers 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Paul  E.  Bishop,  Major,  US  AF 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND‘ADDRESS(ES) 

Air  Force  Institute  of  Technology,  WPAFB  OH  45433-6583 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

AFIT/GE/ENG/90D-05 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION  /AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

A  comprehensive  study  of  the  effects  of  noise  jamming  on  various  pulsed  radar  receiver  types  is  proposed.  First,  a 
literature  review  on  the  sources  of  electrical  noise  and  various  jamming  waveforms  is  presented.  The  power  spectra  for 
three  noise  jamming  waveforms  are  rigorously  derived.  These  cases  are  shot  noise  in  a  parallel-plate  diode,  a  series  of 
pulses  with  random  amplitudes  and  intervals,  and, for  a  series  of  pulses  with  random  amplitudes,  spacing  and  phases.  Next, 
the  probability  density  functions  for  the  output  of  linear  and  quadratic  detectors  are  developed.  For  each  detector  type, 
the  probability  density  functions  for  the  cases  of  noise  only  and  signal  plus  noise  are  derived.  Finally,  the  effects  of  CW 
interference  on  conventional  pulsed  radar  receivers  is  accomplished.  This  analysis  shows  that  the  effect  of  a  pure  CW  tone 
as  an  interfering  signal  is  that  it  increases  the  minimum  power  required  for  detection.  As  the  CW  tone  is  detuned  from  the 
signal  carrier  frequency  the  minimum  power  required  for  detection  increases  until  the  detuning  is  so  great  that  the  effects 
of  the  CW  is  negligible.  This  thesis  lays  the  ground  work  for  a  much  broader  future  study  of  the  effects  of  noise  jamming 
on  various  pulsed  radar  receiver  types. 

14.  SUBJECT  TERMS 

Power  Spectra,  Probability  Density  Functions,  Radar  Jamming,  Noise,  Noise  (Radar),  Noise 
Jamming 

15.  NUMBER  OF  PAGES 

117 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540*01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  bv  ANSI  Std.  Z39-18 
298*102 


